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I.  Summary 

A  single-pulse  shock  tube  has  been  used  to  examine  the  pyrolysis  and  rich  oxidation  of  cy- 
clopentadiene  as  well  as  its  pyrolysis  in  the  presence  of  acetylene,  biacetyl,  and  benzene.  In  addi¬ 
tion,  mixtures  of  benzene  and  hydrogen  have  been  copyrolyzed  at  ratios  of  one-to-three.  These  fuels 
have  been  diluted  in  argon  and  shock  heated  over  the  temperature  range  of  1200  to  2000K  and  at 
total  pressures  of  ten  to  thirteen  atmospheres.  Dwell  times  were  about  500-600  microseconds.  Col¬ 
lected  gas  samples  were  analyzed  using  gas  chromatography  for  hydrogen,  carbon  oxides  and  Ci-  to 
Ci4-hydrocarbons.  Experimental  results  and  preliminary  chemical  kinetic  modeling  lead  to  mech¬ 
anistic  proposals  on  the  decomposition  of  Cs-rings  as  well  as  isomerization  processes  between  C5- 
and  Cc-rings.  In  addition,  data  on  cyclopentadiene  indicates  that  this  fuel  has  a  sooting  tendency 
comparable  to  that  of  benzene.  These  results  indicate  that  Cj-rings  can  play  a  very  significant  role 
in  the  formation  of  soot. 

Simplified  modeling  techniques  of  soot  formation  based  on  the  MAEROS  code  have  been  ex¬ 
tended.  Soot  formation  as  measured  by  Bockhorn  and  coworkers  in  a  premixed  propane/oxygen 
flame  has  been  predicted  although  better  agreement  with  experimental  data  was  obtained  assuming 
that  naphthalene  is  the  inceptor  for  soot.  Previous  studies  had  indicated  that  benzene  formation 
could  be  used  as  a  surrogate  for  the  inception  rate  of  soot  particles.  In  addition,  some  parametric 
studies  on  the  effect  of  variation  in  temperature  and  inception  rates  have  been  performed.  Calcula¬ 
tions  indicate  that  initial  soot  formation  rates  are  strongly  dependent  on  inception  rates,  but  final 
soot  fractions  produced  in  flames  are  much  less  dependent  on  inception  rates  than  had  previously 
been  assumed. 

II.  Introduction 

There  is  increeising  experimental  evidence  that  important  rate-limiting  steps  to  soot  formation 
are  the  production  and  growth  of  aromatic  rings.  Detailed  modeling  and  comparison  to  experimental 
data  has  led  to  a  good  understanding  of  mechanisms  and  rates  for  the  production  of  benzene  (and 
phenyl  radical).  Existing  mechanisms  for  the  formation  of  naphthalene  (plus  related  species)  and 
higher  order  aromatics  are  principally  thermochemical  estimates  and  no  experimental  verification 
of  these  steps  (mechanisms  or  rate  constants)  are  available.  A  major  unknown  is  the  importance  of 
Cs-species  on  ring  growth.  Flame  studies  indicate  that  cyclopentadiene  has  sooting  characteristics 
similar  to  that  of  aromatics,  but  mechanisms  are  lacking  to  explain  the  rapid  conversion  to  Ce- 
rings  which  are  believed  to  dominate  in  ring  growth  processes.  Another  large  uncertainty  in  soot 
formation  mechanisms  is  the  effect  of  oxidation.  Under  many  conditions,  growth  of  polyaromatic 
hydrocarbons  (PAH’s)  may  be  quite  fast  but  is  counterbalanced  by  oxidation.  Therefore,  the  net 
production  rate  of  PAH’s  and  therefore  of  soot  is  only  a  small  fraction  of  the  “gross”  production 
rate  of  PAH’s.  Oxidation  processes  (mechanisms  and  rates)  of  aromatic  hydrocarbons  at  elevated 
temperatures  (>  1200K)  are  unfortunately  very  poorly  known. 

To  investigate  these  important  and  competitive  phenomena,  an  experimental  and  modeling 
program  is  being  performed.  The  experimental  phase  of  the  program  involves  the  use  of  a  single¬ 
pulse  shock  tube  to  thermally  stress  mixtures  of  gases,  which  are  then  quenched,  and  the  reaction 


1 


products  are  analyzed.  This  program  is  focusing  on  detection  of  two-  and  three-ringed  products 
in  order  to  test  mechanisms  proposed  for  the  formation  of  these  species.  In  support  of  this  work 
and  under  corporate  sponsorship,  the  gas  sampling  system  has  been  modified  to  collect  and  detect 
high  molecular  weight  species.  Results  of  these  modifications  were  extremely  successful  and  the 
modified  facilities  were  described  in  last  year’s  report.  To  examine  the  importance  of  C5  rings, 
cyclopentadiene  heis  been  pyrolyzed,  oxidized  and  copyrolyzed  with  acetylene,  biacetyl,  and  benzene. 
Results  on  rich  oxidation  of  benzene  heis  been  investigated  and  were  reported  in  last  years  annual 
report.  Furthermore,  in  order  to  help  obtain  additional  information  on  species  identification,  a  mass 
spectrometer  will  be  added  on  to  the  gas  chromatogragh.  This  added  feature  will  insure  greater 
confidence  in  the  identification  of  high  molecular  weight  intermediates. 

Two  types  of  modeling  are  being  performed  as  part  of  this  contract.  Detailed  chemical  kinetic 
modeling  is  being  used  for  comparison  to  experimental  data  and  to  test  reaction  mechanisms  and 
rate  constants  for  formation  of  multi-ringed  species  and  for  oxidation  processes.  Secondly,  using 
concepts  developed  from  this  work  and  from  the  literature,  simplified  models  are  being  examined 
which  may  be  useful  in  describing  the  formation  and  growth  of  soot.  Preliminary  results  of  calcu¬ 
lations  described  last  year  for  an  ethylene/air  flame  have  been  extended  to  include  calculations  of 
particulates  produced  in  other  laboratory  flames. 

III.  Results 

A.  Experimental  Results  and  Interpretations 

Since  the  end  of  the  previous  AFOSR  contract  in  November  1987,  thirty-eight  separate  mixtures 
have  been  shock  heated  and  products  analyzed  by  gas  chromatography.  The  mixtures  are  listed  in 
Table  I  (for  1988)  and  Table  II  (for  1989).  Six  of  the  mixtures  were  examined  cis  part  of  the  first 
year  of  this  AFOSR  program  and  10  mixtures  were  shock-heated  in  the  second  year  of  the  program; 
the  remainder  were  tested  under  corporate  sponsorship.  For  these  tests,  dwell  times  were  typically 
500  microseconds  and  total  pressures  (balance  argon)  were  9  to  13  atmospheres.  Typically  each 
series  of  experiments  consisted  of  8  to  15  individual  shocks. 

For  each  of  these  series,  gas  chromatographic  data  for  hydrogen,  carbon  oxides,  and  hydro¬ 
carbon  species  from  methane  up  to  C^-species  was  obtained  for  product  gases.  Accuracy  was 
decreased  with  increasing  mass.  Methane  was  measured  to  about  three  percent  accuracy,  benzene 
was  measured  to  about  seven  percent  accuracy  and  the  heaviest  hydrocarbons  had  about  thirty 
percent  uncertainty.  The  complexity  of  the  geis  chromatographic  spectrum  led  to  some  uncertainty 
in  product  identification  of  the  high  molecular  weight  hydrocarbons.  Consequently,  the  analytical 
system  is  being  converted  into  a  gc/ms  system  to  enhance  the  qualitative  analysis  of  this  class 
of  species.  The  mass  spectrometer  will  be  equipped  with  chemical  ionization  as  well  as  electron 
ionization  to  assist  in  peak  identification.  It  was  originally  desired  to  set  up  a  capability  to  measure 
species  concentrations  in  real  time  using  a  time-of-flight  mass  spectrometer  coupled  directly  to  the 
shock  tube.  Unfortunately,  due  to  technical  problems,  this  approach  was  terminated. 

1.  Cyclopentadiene  decomposition  and  ring  formation 

Cyclopentadiene  (CPD)  is  not  sold  commercially  since  it  easily  dimerizes  into  dicyclopentadiene 
(DCP).  Initial  shock  tube  experiments  (corporate  funded)  using  DCP  had  ambiguous  interpreta¬ 
tions,  so  a  facility  was  set  up  to  produce  cyclopentadiene  for  subsequent  shock  tube  experiments. 
DCP  was  decomposed  to  CPD  at  ITO^C,  the  boiling  point  of  DCP.  Subsequently,  CPD  weis  sep¬ 
arated  from  DCP  in  a  distillation  column  and  then  condensed.  Purity  of  the  CPD  was  approx¬ 
imately  96%  which  could  be  further  refined  to  99.6%  by  bulb-to-bulb  distillation  when  sufficient 
CPD  WM  produced.  The  two  principal  impurities  eluted  just  before  CPD  and  are  believed  to  be  C5- 


Table  I 


Series  of  Experiments  Completed  During  the  First  Year  of  AFOSR  Program 


Reactants 

Initial  Concentrations 

Toluene 

1% 

Benzene 

1.1% 

Benzene/oxygen 

1.1/0.22% 

Benzene/oxygen 

0.11/0.022% 

Benzene/oxygen 

1. 1/1.1% 

Dicyclopentadiene 

0.125% 

*Acetylene/ethene 

1.1/2. 2% 

*Acetylene/ethene 

0.125/0.25% 

*Toluene/methanol 

1.0/0.15% 

*  Toluene/hydrogen 

1. 0/3.1% 

*Methylcyclohex3ne 

1% 

*Methanol 

4% 

*Natural  Gas/oxygen 

4/1% 

*Methane/oxygen 

4/1% 

*Cyclopentadiene 

1.4% 

*Cyclopentadiene 

0.1% 

Performed  under  corporate  sponsorship 
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Table  II 


Series  of  Experiments  Completed  During  the  Second  Year  of  AFOSR  Program 


Reactants 

Initial  Concentrations 

CPD"'' /oxygen 

1.4%/0.91% 

CPD/oxygen 

0.1%/0.065% 

CPD/benzene 

1.4%/0.4% 

CPD/benzene 

0.1%/0.03% 

CPD/biacetyl 

1.4%/0.35% 

CPD/biacetyl 

0.25%/0.062% 

CPD/acetylene 

1.4%/ 1.05% 

CPD/acetylene 

0.14%/0.105% 

Benzene /hydrogen 

l.l%/3.3% 

Benzene/hydrogen 

0.335%/!% 

*Methane/oxygen 

4%/0.5% 

*Methane/oxygen 

0.6%/0.075% 

*ToIuene/oxygen 

1%/1% 

*Toluene/oxygen 

0.255%/!% 

*  Heptane 

1% 

*  Heptane 

0.1% 

*Heptane/oxygen 

1.0%/0.5% 

*Heptane/oxygen 

0.1%/0.05% 

■'"cyclopentadiene 

*  performed  under  corporate  sponsorship 
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compounds.  Based  on  retention  times,  the  impurities  have  been  tentatively  identified  as  1-pentene 
and  1,3-pentadiene. 

Including  corporate-funded  experiments,  a  total  of  eight  series  of  experiments  on  CPD  were 
completed.  These  included  pyrolysis  (although  a  small  oxygen  impurity  was  present),  oxidative 
pyrolysis,  and  copyrolyses  with  biacetyl  (a  source  of  methyl  radicals)  and  with  acetylene. 

Interpretation  of  these  data  hzis  been  complex  due  to  a  large  number  of  high  molecular  weight 
species,  the  lack  of  any  literature  mechanisms  for  the  decomposition  pathways,  and  unknowns  re¬ 
garding  the  role  of  Cj-species  in  ring  formation  or  growth.  The  overall  decomposition  of  cyclopen- 
tadiene  for  1.4%  cyclopentadiene  with  various  additives  is  shown  in  Figure  1.  Pyrolysis  exhibits 
the  slowest  rates  of  decomposition,  while  the  copyrolysis  with  biacetyl  exhibits  the  highest  rates 
presumably  because  of  the  high  rates  of  initiation  from  biacetyl^.  Plots  of  the  light  hydrocarbon 
species  as  a  function  of  initial  post-shock  temperatures  are  provided  in  Figures  2  and  3  for  the 
pyrolysis  and  the  oxidative  pyrolysis.  The  distribution  of  products  appears  similar  although  the 
concentrations  produced  during  the  pyrolysis  attain  relatively  lower  levels  (note  the  difference  in 
scales).  Of  particular  note  is  that  the  oxidation  produces  substantially  greater  relative  concen¬ 
trations  of  1,3-butadiene  (1,3-C4H6)  and  vinylacetylene  (C4H4).  It  is  believed  that  these  products 
arise  from  an  oxidation  process  of  the  cyclopentadienyl  radical  similar  to  that  previously  proposed^. 
This  decomposition  route  which  includes  the  formation  of  C4-hydrocarbons  differs  from  a  proposed 
pyrolytic  process  presented  in  the  next  paragraph. 

A  preliminary  mechanism  for  the  pyrolysis  is  presented  in  Table  III.  This  proposed  mechanism 
does  not  include  any  growth  reactions  and  features  decomposition  which  is  preceeded  by  H-atom  ad¬ 
dition  to  CPD.  The  proposed  decomposition  of  the  radical  adduct  to  acetylene  and  the  allyl  radical 
is  perhaps  the  lowest  energy  decomposition  pathway  for  the  Cs-ring.  The  overall  endothermicity  of 
Reaction  3  is  only  about  10  kilocalories  per  mole  although  the  barrier  to  the  reaction  was  assigned 
to  that  for  H-atom  abstraction  from  olefinic  hydrocarbons.  Other  possible  reactions  involve  H-atom 
shifts  within  the  CPD  molecule  followed  by  isomerization  to  a  linear  aliphatic.  Such  steps  are  sim¬ 
ilar  to  those  proposed  for  the  decomposition  of  pyrrole®.  Results  of  calculations  using  this  scheme 
are  shown  in  Figure  4  and  compare  qualitatively  well  with  the  experimental  profiles  of  acetylene 
and  methane  in  Figure  2.  Other  products  not  predicted  by  this  mechanism  but  shown  in  Figures 
2  and  3  may  be  from  the  direct  decomposition/isomerization  of  CPD,  from  decomposition  of  the 
hydrocarbon  impurities,  or  perhaps  from  a  small  fraction  of  oxidation  due  to  an  oxygen  impurity. 
This  preliminary  mechanism  underpredicts  the  amount  of  decomposition  into  products  at  elevated 
temperatures.  Probably  unimolecular  decomposition  or  isomerization  of  CPD  or  the  cyclopentadi¬ 
enyl  radical  occurs  more  rapidly  at  elevated  temperatures.  At  HOOK,  the  model  predicts  that  more 
than  10%  of  the  CPD  is  converted  to  the  cyclopentadienyl  radical  which  is  relatively  unreactive 
according  to  the  preliminary  model.  This  prediction  is  probably  not  accurate.  The  fate  of  the 
cyclopentadienyl  radical  will  have  to  be  investigated  further  although  it  is  reasonable  to  eissume 
that  it  plays  a  significant  role  in  ring  growth  reactions. 

An  important  feature  of  the  data  from  CPD  pyrolysis  is  the  large  amount  of  mass  that  i.s 
converted  to  high  molecular  weight  material.  A  plot  of  the  carbon  contained  in  all  species  with 
molecular  weights  of  128  gram/mole  and  less  is  shown  Figure  5.  This  data  is  compared  with  that 
observed  during  the  pyrolysis  of  other  hydrocarbons.  This  figure  demonstrates  that  cyclopentadiene 
produces  large  quantities  of  high  molecular  weight  material  (a  qualitative  measure  of  soot  produc¬ 
tion)  similar  to  that  of  aromatic  hydrocarbons  and  much  greater  than  observed  for  aliphatic  species. 
This  result  is  not  surprising  considering  the  aromatic  character  of  CPD  and  the  resonantly-enhanced 
stability  of  the  cyclopentadienyl  radical.  This  result  however  leads  to  some  concerns  about  existing 
models  for  soot  production  since  none  of  them  include  the  potential  role  of  Cg-species. 
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Light  Products  from  Pyrolysis  of  1.4%  Cyclopentadiene 
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Table  III 

Preliminary  Mechanism  for  the  Pyrolysis  of  Cyclopentadiene 


Reactions  Considered 

Pre  Exp 

Act  Eng 

1. 

c-CsHe^  H+c-CgHs 

0.600D+16 

81000. 

2. 

H+c-CsHe^  H2+C-C5H5 

O.lOOD+14 

3000. 

3. 

H+c-CgHe'^  C2H2+C3H5 

0.500D+14 

14000. 

4. 

H+c-CjHe'^  C2H2+C3H3+H2 

0.300D+14 

27000. 

5. 

CH2+C2H2''^  C3H3+H 

0.180D+13 

0. 

6. 

CH3  +  C2H2-^  CH3CHCH 

0.620D+12 

7700. 

7. 

H  +  C3H4  CH3CHCH 

0.580D+13 

3100. 

8. 

CHaCHCH^  C3H5 

0.140D+14 

36000. 

9. 

H+ALLENE4->  C3H5 

0.400D+13 

2700. 

10. 

CH3+C-C5H6-^  CH4+C-C5H5 

0.500D+13 

5000. 

11. 

CH3+H--  CH2+H2 

0.724D+15 

15100. 

12. 

C3H3  +  H2’^  C3H4+H 

0.400D+13 

18000. 

13. 

C3H3+C'C5H6<— *  C3H4+C-C5H5 

O.lOOD+14 

8000. 

14. 

CsHs+c-CsHs*-*  AI+C2H2 

0.400D+13 

0. 

15. 

H  +  C2H2  C2H3 

0.550D+13 

2500. 

16. 

C2H4  +  H  C2H3  +H2 

0.692D+15 

14500. 

17. 

CsHs+H'*— *  CH3+C2H3 

0.200D+14 

12000. 

18. 

H+CaHe^  C2H4+CH3 

O.lOOD+13 

5000. 

19. 

CaHe+H^-*  C3H5+H2 

0.501D+14 

3500. 

20. 

C3H5-*  C3H4+H 

0.398D+14 

70000. 

21. 

C3H5  +  H-  CgHs 

0.400D+14 

0. 

Notes:  Units  in  cc,  moles,  cal 

C-C5H5  =  cyclopentadienyl  radical 
c-CjHp  =  cyclopentadiene 
A1  =  benzene 
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1.4%  Cyclopentadiene  Pyrolysis 


Figure  4 
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Mass  Recovery  Following  Pyrolysis  in  a 
Single-Pulse  Shock  Tube 


Figure  5 
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Since  aromatic  species  have  been  linked  to  the  growth  to  soot,  it  was  of  interest  to  compare 
benzene  production  from  cyclopentadiene  to  that  produced  from  other  hydrocarbons.  Benzene 
production  has  been  previously  shown  to  correlate  well  with  soot  production^  ®.  Shown  in  Figure 
6  is  the  benzene  production  from  cyclopentadiene  as  well  as  that  from  a  few  other  fuels.  As  ran 
be  seen  in  this  figure,  the  level  of  benzene  production  can  be  considered  intermediate  relative  to 
that  produced  from  other  fuels,  yet  the  mass  deficit  apparent  from  Figure  5  is  much  larger  than 
the  other  fuels  included  in  this  figure.  Thus,  benzene  production  from  cyclopentadiene  does  not 
correlate  with  soot  production  and  suggests  that  benzene  formation  is  skipped  in  the  production 
of  high  molecular  weight  species.  These  results  regarding  ring  growth  from  CPD  are  significant 
since  they  demonstrate  that  Cs-rings  (as  well  as  polyaromatic  hydrocarbons  containing  Cj-rings) 
are  much  more  important  then  present  models  describing  soot  formation  assume. 

Rapid  interchange  between  C5  and  Ce-rings  may  also  occur.  To  test  this  hypothesis,  a  series 
of  pyrolysis  experiments  were  performed  with  additives.  One  of  the  additives,  biacetyl  ((CH3)2CO) 
was  selected  since  it  provides  a  facile  source  of  methyl  radicals  and  easy  modeling  techniques  are 
available  for  estimating  the  methyl  radical  concentrations*.  Methyl  radicals  in  turn  may  add  to 
the  cyclopentadienyl  radical  to  form  methylcyclopentadiene.  After  loss  of  an  H-atom,  the  resulting 
radical  may  isomerize  to  cyclohexadienyl  lose  another  H-atom  and  eventually  form  benzene.  This 
process  is  similar  to  the  reverse  of  the  mechanism  for  benzene  decomposition  recently  proposed  by 
Ritter,  et  al.^.  A  second  mechanism  for  ring  interchange  wets  investigated  by  adding  acetylene  to  the 
pyrolysis  of  cyclopentadiene.  We  speculated  that  acetylene  add  to  cyclopentadienyl  radical  and  the 
resultant  adduct  undergoes  a  ring  enlargement  reaction  to  form  the  benzyl  radical,  which  becon.es 
toluene  after  H-atom  addition.  This  process  is  the  reverse  of  a  proposed  step  for  the  decomposition 
of  the  benzyl  radical^,  although  many  different  proposals  have  been  made  (see  for  example  Ref.  8). 

Benzene  production  from  pure  cyclopentadiene  is  compared  to  benzene  production  during  copy¬ 
rolysis  of  cyclopentadiene  and  biacetyl  in  Figure  7.  Also  included  on  this  plot  is  data  on  the  two 
isomers  of  the  speculated  intermediate,  methylcyclopentadiene.  Biacetyl  clearly  and  dramatically 
increases  the  concentration  of  the  methylcyclopentadienes  as  well  as  significantly  increase  the  con¬ 
centration  of  benzene.  In  Figure  8,  toluene  produced  from  cyclopentadiene  pyrolsis  and  from  a 
mixture  of  CPD  and  acetylene  is  shown.  Again  large  increases  in  the  conversion  to  aromatics  was 
observed  by  the  presence  of  an  additive.  Together  these  results  not  only  provide  evidence  support¬ 
ing  recent  proposals  for  ring  fracture,  but  also  confirm  speculation  of  rapid  interchange  mechanisms 
between  C5  and  Cc, -rings. 

2.  Explanation  for  Methane  Formation  During  Benzene  Pyrolysis 

In  many  previous  experiments  on  benzene  pyrolysis,  a  significant  amount  of  methane  has  been 
produced  (on  the  order  of  0.5%  of  the  carbon)®.  The  magnitude  of  this  methane  production  is  large 
considering  existing  models  for  benzene  production  predict  either  no  or  negligible  conversion  to 
methane.  The  observed  formation  of  methane  has  remained  unexplained  for  several  years.  The  work 
(experimental  and  modeling)  described  above  on  cyclopentadiene  pyrolysis  offers  a  very  reasonable 
explanation  for  this  methane  formation.  The  general  reaction  scheme  is  as  follow?,  although  nnl.\  a 
portion  of  the  benzene  («  5%)  decomposes  in  this  manner 

H  +  CeHc,  -  CcHt 
CoHt  --  MeCPD 
MeCPD-bH^-CHg  +  CsHs 
CH3  +  CeHe  -  CH4  +  CeHs 
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Benzene  Production  During  Pyrolysis  in  a 
Sinde-Pulse  Shock  Tube 


Figure  6 
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Formation  of  Benzene 
and  Methylcyclopentadienes 


Figure  7 
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Formation  of  Toluene 


Figure  8 
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Subsequent  decomposition  of  the  cyclopentadienyl  radical  (C5H5)  can  be  expected  to  form  some 
additional  methane.  At  the  temperatures  of  benzene  pyrolysis  (1600-1700K),  this  decomposition 
should  be  fast. 

S.  Addition  of  Hydrogen  to  Benzene  Pyrolysis 

Profiles  of  end-product  species  produced  from  the  copyrolysis  of  1.1%  benzene/3.3%  hydrogen 
were  compared  with  previously  obtained  data  on  the  pyrolysis  of  1.35%  benzene.  Of  particular 
note  was  the  fact  that  hydrogen  addition  suppressed  both  the  decomposition  of  benzene  at  lower 
temperatures  and  the  formation  of  high  molecular  weight  species.  Both  of  these  effects  can  be 
traced  to  the  reaction 


C6H5  +  H2  --CeHo  +  H 

which  is  pushed  to  the  right  with  excess  hydrogen  and  thus  suppresses  the  phenyl  radical  concentra¬ 
tion.  Phenyl  radical  in  turn  is  the  controlling  intermediate  species  which  leads  to  ring  decomposition 
as  well  as  to  ring  growth.  Plots  of  some  polyaromatic  hydrocarbons  are  shown  in  Figure  9  for  the 
benzene/hydrogen  mixture  and  can  be  compared  to  that  produced  from  the  pyrolysis  in  Figure  10. 
Hydrogen  addition  leads  to  approximately  a  30%  reduction  in  biphenyl  and  about  a  factor  of  two 
reduction  in  acenaphthylene.  Phenanthrene  production  is  essentially  fully  suppressed. 

A  simplified  model  was  constructed  for  benzene  pyrolysis  based  on  reactions  and  rate  constants 
previously  published^®,  although  the  rate  constant  by  Fahr  and  Stein''  for  the  reaction 

CeHs  +C6H6  Ci2H,o+  H 

was  used.  Predictions  of  the  biphenyl  production  are  shown  in  Figure  11.  Calculations  are  shown 
for  both  the  pyrolysis  and  the  copyrolysis  with  hydrogen.  Qualitatively,  the  predictions  are  very 
similar  to  the  experimental  data.  Perhaps  the  most  noticeable  differences  is  the  location  of  the 
peak  which  is  predicted  to  be  about  lOOK  higher  than  the  peak  in  experimental  data.  In  addition 
the  experiment  indicates  substantial  quantities  of  biphenyl  formation  at  low  temperatures.  Using 
a  preliminary  mechanism  for  the  oxidative  pyrolysis  of  benzene,  both  of  these  differences  could  be 
explained  by  small  concentrations  (a  few  hundred  ppm)  of  oxygen.  The  absolute  magnitude  of 
biphenyl  is,  however,  not  dramatically  effected. 

In  modeling  work  on  toluene,  it  was  found  that  substantial  quantities  of  bibenzyl  were  formed 
in  the  rarefraction  wave  in  the  shock  tube  prior  to  collection  and  analysis  of  the  gas  sample.  To 
test  whether  a  similar  phenomenon  was  occurring  with  biphenyl  in  the  benzene  system,  calculated 
profiles  of  biphenyl  as  function  of  time  were  plotted.  These  curves  are  shown  in  Figure  12.  At 
low  temperatures,  the  production  of  biphenyl  is  dramatically  reduced  once  the  rarefraction  wave 
arrives  (at  600  microseconds).  This  occurs  even  at  elevated  temperatures  when  then  the  biphenyl 
roncpntration  increases  to  a  maximum  and  then  fall'  befnrp  the  arrival  of  the  cooling  wave.  Thus 
the  rarefraction  wave  appears  to  quench  the  production  and  decay  nf  biphenyl  so  that  the  measnrnd 
concentrations  accurately  reflect  the  actual  concentrations  at  the  end  of  the  high  temperature 
pyrolysis  period. 

B  Soot  Aerosol  Dynamics  Simulations 

In  the  last  annual  report,  we  presented  the  results  of  a  simplified  aerosol  dynamics  model 
for  soot  growth.  With  the  aim  of  modelling  the  growth  of  soot  spheroids,  we  adapted  the  most 
recent  version  of  a  well-known  aerosol  dynamics  code,  MACROS,  which  is  a  fixed-  grid,  sectional 
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model  including  nucleation,  surface  growth,  and  coagulation.  The  lowest  size  class  was  assumed  to 
approximate  that  of  a  benzene  molecule,  with  an  effective  nucleation/inception  rate  given  by  the 
production  rate  of  the  smallest  aromatic  molecules,  typically  taken  to  be  benzene.  Surface  growth 
was  assumed  to  be  due  to  acetylene  deposition  alone,  at  a  rate  given  by  tlie  Harris- Weiner  growth 
rate  expression'^.  Provision  for  oxidative  mass  removal  at  a  rate  given  by  the  Nagle  and  Strickland- 
Constable  expression was  also  made  in  the  model.  The  most  important  approximation  inherent 
in  this  approach  is  the  extrapolation  of  particle  dynamic  processes  and  rates  into  the  pre-particle 
regime.  This  is  a  simplified  picture  of  soot  growth  that  glosses  over  many  details  of  the  inception 
process;  the  only  medium  parameters  in  addition  to  temperature  that  need  to  be  specified  are 
the  inception  (benzene  formation)  rate  and  the  concentration  of  the  condensible  acetylene  vapor. 
The  first  comparison  of  the  model  with  the  experimental  volume  fraction  data  of  Harris  and  co- 
workers'^''*,  taken  in  a  premixed,  atmospheric  pressure,  ethylene-oxygen  flame,  was  very  good. 
Both  the  magnitude  and  time  dependence  of  the  soot  growth  were  modelled  well,  and  there  wa 
no  need  to  allow  for  oxidation,  a  result  that  was  consistent  with  the  stoichiometry.  Given  the 
simplifications  made,  we  recognized  the  possibility  that  the  agreement  obtained  was  fortuitous,  and 
initiated  efforts  to  model  other  well-characterized  flames.  One  set  of  calculations  in  progress  is  to 
see  if  the  model  reproduces  the  stoichiometric  variation  of  soot  production  in  the  Refs.  12,14  flames; 
this  work  involves  using  CHEMKIN  to  calculate  the  acetylene  concentration  and  the  benzene  source 
rate,  and  is  in  progress  at  this  time. 


The  only  other  well-characterized  data  suitable  for  modelling  purposes  have  been  reported 
by  Bockhorn,  et.al.""  in  low-pressure,  premixed  flames.  Soot  volume  fraction  profiles  and  all  the 
other  input  data  needed  to  perform  the  modelling  calculations  are  available.  The  volume  fraction 
data  extend  over  a  temperature  range  of  about  450  K  for  the  propane-oxygen  flame  we  initially 
concentrated  on.  The  Maeros2  program  which  we  use  employs  a  linear  interpolation  scheme  to 
determine  coagulation  and  condensation  coefficients,  these  having  been  calculated  exactly  at  four 
combinations  involving  two  input  temperatures  and  pressures.  We  found  that  the  450  K  interval 
was  too  large  for  accurate  linear  interpolation,  and  thus  had  to  confine  attention  to  the  first  ten 
milliseconds.  Work  has  begun  on  modifications  to  the  program  to  permit  interpolation  between 
four  input  temperatures,  which  should  suffice  for  the  larger  temperature  variations.  In  any  event, 
considerable  growth  of  the  soot  volume  fraction  has  occurred  within  the  first  ten  milliseconds,  and 
it  is  not  felt  that  any  incorrect  conclusions  will  result  from  this  restriction  on  the  calculations.  The 
Harris- Weiner  surface  growth  rate  was  modified  to  include  a  31.8  kcal  activation  energy'^',  and  the 
oxidation  rate  has  the  temperature  dependences  given  in  Ref.  13.  Calculations  performed  so  far 
have  been  limited  to  the  0.15  atmosphere  propane-oxygen  flame  of  Ref.  15,  Figure  13  shows  the 
comparison  of  the  calculated  and  experimental  soot  volume  fractions,  assuming  that  the  benzene 
source  rate  of  Ref.  15  is  the  appropriate  inception  rate.  With  no  oxidation,  the  calculated  volume 
fractions  are  too  large,  but  there  is  considerable  sensitivity  to  the  oxygen  concentration,  as  seen. 
For  an  oxygen  concentration  of  about  1%,  the  agreement  is  reasonable  for  longer  times.  This  flame 
is  somewhat  leaner  than  the  Harris  flame  discussed  e?ir!ipr.  "^o  the  need  to  incorporate  oxidation  is 
at  least  qualitatively  reasonable.  Some  indication  of  (he  predicted  tomperaturf  =;ensi(ivi(y  is  shown 
in  Figure  14.  A  reasonable  temperature  variation  of  ±  100  K  does  not  appear  to  explaiti  the  theory- 
experiment  discrepancy.  The  foregoing  calculations  assumed  a  coalescence  sticking  probability  of 
unitv.  Choosing  a  value  of  2.2,  as  recommended  by  Harris  and  Kennedy''  improves  the  agreement 
some  '-hat,  but  does  not  change  the  basic  conclusions.  A  slightly  smaller  oxygen  concentration  is 
required  for  agreement,  naphthalene  production  rates  are  also  available  in  Ref.  16.  and  we  have 
found  that  taking  naphthalene  as  the  inception  source  results  in  much  better  agreement,  as  shown 
in  Figure  15.  The  naphthalene  source  rate  is  about  one-sixtieth  that  of  benzene.  A  Van  der  Waals 
sticking  probability  of  2.2  has  been  employed  in  this  calculation.  Our  sensitivity  calculations  show 
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Simulation  of  Bockhorn  Flame 


Figure  13 
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Simulation  of  Bockhorn  Flame 


Figure  15 
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that,  compared  to  the  benzene-based  simulations,  smaller  temperature  uncertainties  and  oxygen 
concentrations  would  need  to  be  invoked  to  resolve  the  small  remaining  discrepancies.  Picking 
naphthalene  as  the  inception  size  class  might  seem  to  be  a  better  choice  than  benzene,  inasmuch 
as  it  is  closer  to  the  particle  regime,  but  until  more  information  is  available  about  the  oxygen 
concentrations  in  these  flames,  no  quantitative  conclusions  can  be  drawn.  The  choice  of  benzene 
or  naphthalene  a.s  the  inception  source  mainly  affects  the  calculated  soot  loadings  at  early  times. 
Asymptotically,  at  long  times,  calculated  differences  are  generally  much  less  sensitive  to  inception 
rate.  Because  of  the  previously  discussed  temperature  interpolation  problem,  however,  we  can’t  yet 
present  the  accurate  long  time  behavior  of  the  Ref.  15  simulations.  Some  indication  of  the  sensitivity 
is  presented,  however,  in  Figure  16,  where  for  a  constant  temperature  of  1650  K  the  effect  of  reducing 
the  nominal  (benzene)  inception  rate  by  a  factor  of  100  is  shown.  At  100  milliseconds,  there  is  only 
about  a  factor  of  two  difference  between  the  two  predictions. 

In  addition  to  improving  the  temperature  range  capability  of  Maeros2,  a  need  was  also  identified 
to  improve  the  numerical  integration  scheme  when  there  is  net  soot  loss  due  to  oxidation;  the  present 
scheme  can  become  unstable  when  oxidation  dominates  acetylene  deposition. 

IV.  List  of  Publications 

A  reply  (to  a  comment)  entitled  “On  Impurity  Effects  in  Acetylene  Pyrolysis”  by  M.  D.  Colket,  II. 
B.  Palmer  and  D.  J.  Seery  has  been  accepted  for  publication  in  Combustion  and  Flame.  A  copy  of 
this  article  is  provided  in  Appendix  B. 

An  article  entitled  “Shock  Tube  Pyrolysis  of  Pyridine”  by  J.  C.  Mackie,  M.  B.  Colket,  and  P.F. 
Nelson  has  been  accepted  for  publication  by  the  Journal  of  Physical  Chemistry.  The  final  manuscript 
is  included  in  Appendix  C.  The  research  was  performed  under  corporate  sponsorship. 

A  manuscript  entitled  “Shock  Tube  Pyrolysis  of  Pyrrole  and  Kinetic  Modeling”  by  J.  C.  Mackie.  M. 
B.  Colket.  P.  F.  Nelson  and  M.  Esler  has  been  submitted  to  the  International  Journal  of  Chemical 
Kinetics  for  review.  A  copy  of  the  original  manuscript  is  in  Appendix  D.  The  research  was  performed 
under  joint  sponsorship  from  UTRC  and  the  University  of  Sydney. 

An  article  entitled  “The  High  Temperature  Pyrolysis  of  Cyclopentadiene”  by  M.  B.  Colket  will  be 
written  in  1990-91  and  submitted  to  Combustion  and  Flame. 

A  manuscript  entitled  “A  Simplified  Model  for  the  Production  of  Soot  in  Premixed  Flames”  by 
R  .1.  Hall  and  M.  B.  Colket  will  be  written  in  1990-91  and  submitted  to  Combustion  Science  and 
Technology. 

V.  Meeting  Interactions  and  Presentations 

1.  A  round  table  discussion  on  “Current  Problems  in  Soot  Formation  During  Combustion, 
Especially  the  Mechanism  of  Soot  Formation”  was  held  in  Gottingen,  West  Germany  on 
March  29-30,  1989.  The  meeting  was  organized  by  Professor  H.  Gg.  Wagner  and  was 
attended  by  about  twenty  engineers/scientists  currently  involved  with  understandinp  soot 
formation  phenomena.  Financial  support  was  provided  by  the  Commission  for  Condensa¬ 
tion  Phenomena  of  the  Academy  of  Sciences  in  Gottingen.  West  Germany. 

2.  Department  of  Energy  -  Office  of  Basic  Energy  Sciences  Combustion  Research  Meeting 
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3  A  paper  entitled  “Simplified  Models  for  the  Production  of  Soot  in  a  Premixed  Flame” 
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by  R.  J.  Hall  and  M.  B.  Colket  was  presented  at  the  Eastern  Section  of  the  Combustion 
Institute  on  October  30^*’  -  November  1*‘,  1989.  The  meeting  was  held  in  Albany,  New 
York.  A  copy  of  this  paper  is  provided  in  Appendix  A. 

4.  M  B.  Colket  presented  an  invited  lecture  entitled  “Progress  Towards  Understanding  Soot 
Formation  and  Development  of  Global  Models”  to  the  Diesel  Cooperative  Meeting  held  at 
United  Technologies  Research  Center  on  May  17-18,  1990. 
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Soot  inception  and  growth  kinetics  have  been  studied  using  two  types  of  modeling.  In  studies  of 
inception,  detailed  kinetic  modeling  has  been  performed  to  test  reaction  mechanisms  and  calculate 
rate  constants  for  formation  of  multiring  aromatic  species.  One  such  process  involves  the  stepwise 
acetylene  addition  mechanism 


A,  -(-  Hj  ^  Ai  +  H 

(1) 

Ai  -l-  C2H2  AiC2H  -b  H 

(2) 

AiC2H-l-H^  AiC2H-(-H2 

(3) 

AiC2H  -b  C2H2  — ‘  Ai+1 

(4) 

where  Aj  is  a  closed  ring  aromatic  species  Ai  =  benzene,  Aj  =  naphthalene,  Ai  =  phenyl,  etc.) 

This  sequence  reduces  to  the  mechanism  found  by  Frenklach,  et  al*  to  dominate  for  growth 
from  one  to  two  aromatic  rings.  A  similar  simplified  growth  process  has  also  been  suggested  by 
Frenklach*.  Rate  constants  used  in  the  present  study  are  listed  in  Table  I.  k_i,  kj  were  obtained 
from  References  3  and  4  for  the  case  of  Ai=benzene.  ks  and  k4  were  estimated  by  analogy  to 
k_i  and  kj,  respectively.  The  remaining  rate  constants  were  calculated  from  equilibrium.  Using 
steady-state  approximations  for  AiCjH  and  AiCjH  analytical  solutions  for  the  concentration  of 
Ai+i  as  a  function  of  time  can  be  obtained.  Based  on  these  results^,  characteristic  times  for  the 
production  of  Ai  can  be  determined.  For  the  experimental  conditions  of  Harris  and  co-workers® ■*, 
the  production  of  a  1000  amu  particle  (~20  aromatic  rings)  would  take  10  to  20  milliseconds,  a  time 
considered  long  for  production  of  an  incipient  particle.  We  come  to  the  conclusion  that  either  (1) 
rate  constants  selected  in  Table  1  are  significantly  in  error,  (2)  the  mechanism  is  not  the  dominant 
reaction  path,  or  (3)  the  simplified  reaction  sequence  does  not  describe  adequately  a  more  complex 
path  for  PAH  growth.  The  steady-  state  approximation  has  been  compared  to  standard  integration 
techniques  and  shown  to  be  a  reasonably  simple  assumption.  We  have  also  considered  variation 
of  rate  constants  due  to  size  (collision  cross-section,  number  of  active  sites,  etc.)  and  although 
some  acceleration  in  rate  is  observed,  the  rate  enhancement  appears  insufficient  to  describe  soot 
production  within  a  few  milliseconds  as  experimentally  observed. 

The  possible  role  of  agglomerative  reactions  of  the  form 

Ai  +  Aj  — >  (Ai)(Aj)  =  Aij  (large  i,j) 
and  of  the  reaction  sequence 

Ai  +  Ai  ^  AiAi  +  H 
AiAi  +  H  ^  AiAi  -t-  Hj 
AiAi  +  CiHj^  AsH-t-H 
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has  also  been  considered.  This  series  of  reactions  is  similar  to  the  sequence  found  by  Frenklach  et 
al®,  to  accelerate  soot  formation  during  benzene  pyrolysis. 

Soot  primary  spheroid  growth  has  also  been  modeled  as  an  aerosol  dynt  nics  problem  using 
the  MAEROS  code®.  The  size  range  of  spherical  particles  is  divided  into  discrete  classes,  and 
an  equation  of  motion  for  each  size  class  is  solved  with  terms  for  nucleation  (inception),  surface 
growth,  oxidation,  and  coalescence.  Using  an  experimentally  determined  value  for  the  benzene 
formation  rate®,  calculations  for  hydrogen  atom  concentrations^,  choosing  the  smallest  size  class  to 
be  that  of  the  benzene  molecule,  and  extrapolating  the  Harris- Wiener  rate  expression*®  for  surface 
growth  by  acetylene  deposition  to  this  size  range,  we  obtained  very  good  theoretical  agreement  with 
soot  volume  fractions  measured  in  a  pre-mixed  flame**.  The  uncertainties  in  the  calculation  will 
be  discussed,  and  the  sensitivity  to  such  quantities  as  the  coalescence  sticking  probabilities,  the 
number  of  size  intervals,  and  rate  constant  for  surface  growth  will  be  discussed. 

Figure  1  shows  the  calculated  evolution  of  the  normalized  soot  spheroid  size  distributions 
with  time  for  the  premixed  flame  conditions  of  Ref.  6  and  7.  With  simultaneous  nucleation, 
surface  growth,  and  coalescence  (assumed  to  occur  here  with  unit  sticking  probability),  a  “self¬ 
preserving”  form  is  attained  after  about  four  milliseconds  in  this  atmospheric  pressure  simulation. 
Figure  2  shows  the  calculated  evolution  of  the  soot  volume  fraction,  using  the  coalescence  sticking 
probability  as  a  parameter.  For  slow  coagulation,  the  surface  area  of  the  spheroids  remains  high, 
surface  growth  is  correspondingly  rapid,  and  the  computed  volume  fraction  large.  For  sticking 
probabilities  in  the  range  one  to  2.2,  the  agreement  of  theory  with  the  experimental  measurements 
of  Ref.  11  (interpolated  estimate  for  c/o=0.92)  is  very  good. 


Table  I 


Rate  Constants 

for  Growth  Reactions 

Reaction 

log  (A*for) 

Efor/R 

log  (Arev) 

Erev/R 

1 

12.7 

5400 

14.4 

8052 

2 

13.6 

5083 

15.4 

5885 

3 

14.4 

8052 

12.7 

5400 

4 

13.6 

5083 

*  cm® /mole /sec 
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1. 


Calculated  evolution  of  soot  spheroid  size  distribution  in  simulation  of 


Ref.  6  flame. 


2.  Comparison  of  calculated  ( - )  and  meaisurcd  soot  volume  fraction 

profiles  for  premixed  flame  conditions  of  Ref.  11. 
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ON  IMPURITY  EFFECTS  IN  ACETYLENE  PYROLYSIS 


by 

M.  B.  Colket,  III.,  H.  B.  Palmer*  and  D.  J.  Seery 
United  Technologies  Research  Center 
Ecist  Hartford,  CT  06108 

*  Pennsylvania  State  University- 
University  Park,  PA  16802 


In  our  recent  paper'  (CSP),  we  concluded  that  a  radical-chain  mechanism  of  acetylene  pyrolysis 
must  be  appreciable  and  possibly  dominant  and  argued  that  impurities,  for  example  acetone,  may 
contribute  significantly  to  chain  initiation.  To  support  our  conclusions,  a  model  including  acetone 
chemistry  was  developed  but  it  did  not  invoke  vinylidene  reactions  or  other  non-chain  processes 
despite  reasonable  arguments^  in  support  of  such  steps.  Nevertheless,  we  also  recognized,  and 
stated  in  the  paper,  that  the  vinylidene  mechanism  could  not  be  disproven  and  may  in  fact  play  a 
role. 


This  analysis  apparently  prompted  Duran,  Am.orebieta,  and  Colussi^  (DACl)  to  perform  ad¬ 
ditive  experiments.  They  pyrolyzed  acetone/acetylene  mixtures  at  ratios  of  about  0.2  to  2,  selected 
a  decomposition  rate  expression  and  determined  associated  rate  constants.  Colussi^  extrapolated 
this  rate  expression  down  to  concentration  ratios  of  0.001  and  concluded  that  the  impurity  effect 
can  be  neglected  since  it  can  account  for  only  2  to  14%  of  the  decomposed  acetylene  at  9J0K.  In 
contrast  to  the  suggestion  of  Colussi,  we  do  not  consider  these  results  to  be  in  conflict  with  our 
conclusions. 

First  of  all.  the  uncertainty  in  extrapolating  a  rate  expression  derived  from  experiments  with 
acetone/  acetylene  ratios  of  0.2  down  to  ratios  of  only  0.001  may  be  significant.  Our  calculations  sug¬ 
gest  that  at  high  acetone/acetylene  ratios,  acetone  is  rapidly  depleted  by  radical  (H-atom,  methyl, 
or  vinyl)  attack  and  therefore  is  no  longer  available  to  initiate  the  reaction.  Although  absolute  rate 
constants  for  R  +  acetylene  may  in  fact  be  higher  than  R  +  acetone,  the  former  is  a  reversible 
process  and  the  latter  is  not,  because  of  the  rapid  decomposition  of  the  product  acetonyl  radical 
from  the  latter  reaction.  As  radical  concentrations  increase  (i.e.,  as  the  initial  acetone  concentration 
increases),  the  rate  of  R  -I-  acetone  relative  to  the  rate  of  R  -f  acetylene  increeises  dramatically. 
This  result  differs  from  the  conclusion  reached  in  Appendix  I  of  DACl.  If  acetone  is  rapidly  depleted 
by  radical  attack,  then  radical  initiation  rates  will  be  lower  and  overall  decomposition  rates  will  be 
lower  than  if  acetone  were  not  removed  by  radicals.  Rate  expressions  derived  from  high  concen¬ 
tration  experiments  may  then  underestimate  decomposition  when  extrapolated  to  very  low  acetone 
levels.  Therefore,  we  believe  that  the  experiments  of  DACl  (which  suggest  acetone  could  contribute 
about  10%  to  the  overall  decomposition  of  acetylene)  underestimate  the  role  of  low  concentrations 
of  acetone. 

A  second  reason  for  our  view  of  the  experimental  results  of  DACl  is  that  the  analysis  described 
in  our  paper  suggests  that  acetone  could  play  a  relatively  more  important  role  in  radical  initiation 
at  higher  temperatures.  As  stated  in  the  conclusions  of  our  paper,  900K  is  approximately  the  lower 
limit  for  which  acetone  was  predicted  to  play  a  significant  role.  Our  lower  limit  is  very  close  to  the 
temperature  of  the  DACl  study,  910K. 


1 


The  comment  by  Colussi  implies  that  trace  acetone  concentration  levels  cannot  appreciably  ac¬ 
celerate  overall  rates  of  acetylene  pyrolysis.  This  statement  is  not  apparent  from  the  pre-publication 
manuscript  (DACl),  since  acetone  and  acetylene  were  copyrolyzed  at  comparable  concentrations. 
Typical  concentrations  were  51.4  torr  acetylene  and  100  torr  acetone.  Low-level  impurity  experi¬ 
ments  apparently  were  not  performed. 

We  have  performed  some  calculations  bcised  on  the  conditions  cited  by  Duran,  Amorebiela, 
and  Colussi^  (DAC2)  for  pyrolysis  of  neat  acetylene.  We  have  slightly  modified  our  published 
mechanism  to  account  for  the  pressure  dependence.  The  modifications  include  selection  of  k/koo  = 
0.5  for  the  following  reactions. 


H-t-CiHj -- C2H3  (2) 

CH3COCH3  CH3  +  CH3CO  (5) 

In  addition,  the  high  pressure  rate  constant  recommended  by  Benson  and  O’Neal®  is  used  for 
Reaction  5.  At  910K  and  100  torr  acetylene,  we  calculate  a  steady-state  decomposition  rate  (after 
a  transient  less  than  two  seconds)  to  be  2.2x10“®  moles/cm®/sec  (2.2x10”®  M  s”*)  and  a  rate 
constant  of  7.1  x  10®  cm® /moles/sec  (7.1  M”*  s”').  These  values  are  very  close  to  the  values 
cited  by  DACl  and  DAC2.  The  excellent  agreement  is  probably  fortuitous,  but  these  calculations 
were  performed  using  a  chain  mechanism,  an  impurity  of  0.1  torr  of  acetone,  and  did  not  include 
contributions  of  vinylidene.  Detailed  product  information  is  not  available  from  these  publications 
for  comparison. 

It  should  be  emphasized  that  the  thermochemistry  for  several  radical  intermediates  potentially 
important  to  acetylene  decomposition  cannot  be  considered  to  be  well  known.  Consequently  it  is 
difficult  to  draw  definitive  conclusions  baised  on  such  arguments.  The  uncertainty  in  the  heat  of 
formation  of  C4H3  isomers  was  cited  as  an  example  in  CSP,  but  uncertainties  in  the  entropy  also 
exist.  Values  vary  depending  on  the  estimation  (calculation)  technique  as  well  as  the  individual 
performing  the  evaluation.  Most  evaluations  are  unpublished,  but  examples  of  some  reported  values 
of  AH®2<)8>  Sjgg,  and  Cp(T)  for  the  t-C4H3  radical  are  reproduced  in  Table  1.  Also  included  are 
calculations  of  the  rate  constant  for  the  suggested  reaction 

C2H2  +  C2H2  <-+  i  —  C4H3  -l-  H  (f^) 

After  DACl,  kia’s  were  calculated  assuming  the  reverse,  recombination  rate  constant  is  5  x  10^® 
cm®/mole/sec.  DACl  estimated  that  Reaction  la  had  to  have  rate  of  at  leeist  85  cm®/moIe/sec 
at  910K  in  order  for  it  to  contribute  to  acetylene  pyrolysis.  As  can  be  seen  from  the  different 
values  listed  in  Table  I,  any  conclusion  can  be  drawn,  although  most  recent  estimates  using  group 
additivity  techniques  support  the  conclusion  of  DACl  regarding  the  importance  of  Reaction  la. 

Lacking  sufficiently  accurate  thermochemical  data,  one  must  consider  experimental  evidence. 
It  is  our  own  interpretation*  that  available  data  supports  contribution  from  a  chain  mechanism. 
Radical  initiation,  if  not  from  Reaction  la  or  similar  step,  could  arise  from  impurities  or  walls. 
Although  we  have  argued  that  acetone  plays  an  important  role,  acetone  is  not  always  present. 
Acetylene  made  from  water  and  calcium  carbide,  for  example,  will  not  contain  acetone,  although 
other  impurities  (possibly  inorganic)  could  be  present.  If  the  impurity  has  a  relatively  weak  bond, 
large  concentrations  are  not  required. 
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The  difficulty  of  accurately  making  measurements  of  acetylene  pyrolysis  near  lOOOK  should  not 
be  underestimated.  In  addition  to  being  susceptible  to  radical  sources  from  impurities  (see  CSP), 
the  pyrolysis  is  extremely  exothermic.  Formation  of  vinylacetylene,  benzene,  or  solid  carbon  plus 
hydrogen  is  19,  48,  and  54  kilocalories  exothermic  per  mole  of  acetylene  decomposed,  respectively. 
Thus,  nearly  ninety  percent  of  the  potential  energy  release  arises  just  by  forming  benzene  and 
helps  to  explain  the  explosive  behavior  of  neat  acetylene.  Assuming  that  only  2%  of  the  acetylene 
forms  benzene  (an  important  low  temperature  product)  and  taking  Cp(C2H2)  =16  cal/mole/K,  the 
temperature  rise  is  60K  for  undiluted  samples  of  acetylene.  For  a  reaction  with  E=40  kcal/mole,  the 
rate  constant  would  rise  by  a  factor  of  2.5.  For  larger  extents  of  reaction,  this  problem  will  be  more 
severe.  Formation  of  vinylacetylene,  benzene  and  other  higher  molecular  weight  species  will  lead  to 
a  reduction  in  the  number  of  moles  and  to  a  pressure  reduction  (for  constant  volume  reactors)  as  the 
reaction  progresses.  The  pressure  cannot  be  directly  related  to  the  rate  of  reaction  unless  detailed 
information  on  product  distribution  is  also  known.  An  additional  complication  in  interpreting 
pressure  traces  arises  from  the  temperature  rise  due  to  the  heat  of  reaction.  The  temperature  rise 
would  lead  to  a  pressure  increase  (ideal  gas  law)  which  partially  offsets  the  pressure  decay  due  to 
reaction.  For  flow  reactors  in  wh'ch  constant  pressure  can  be  maintained,  the  above  effects  lead  to 
changes  in  velocity  which  complicate  determination  of  reaction  time.  On-line  diagnostics  such  eis 
mass  spectrometers  may  be  pressure  sensitive  and  so  calibrations  and  experiments  must  be  carefully 
interpreted.  In  some  reactor  designs,  energy  releeise  from  the  reaction  may  be  absorbed  by  the  walls, 
but  in  such  situations,  the  rate  of  wall  collisions  must  be  high  and  there  may  be  surface  contributions 
to  the  chemistry.  Pretreatment  of  vessels  at  elevated  temperatures  may  simply  ‘stabilize’  the  source 
of  radical  production/termination  at  the  walls.  The  difficulty  in  accurately  measuring  the  rate  of 
pyrolysis  of  neat  acetylene  is  dramatically  illustrated  by  the  scatter  (factor  of  two  to  three)  obtained 
in  the  work  of  DAC2  (Fig.  1)  but  is  significantly  less  (s»  20%)  when  a  stable  source  of  radicals 
(acetone  or  neopentane)  is  added  (see  Figs.  1-3  in  DACl).  The  scatter  in  the  ‘neat’  acetylene 
experiments  may  be  due  to  run-to-run  variations  in  impurities  or  wall  effects. 

As  we  consider  the  varying  evidence  and  conclusions  from  previous  work,  the  need  for  careful 
documentation  in  future  experiments  becomes  very  evident.  Levels  of  impurities  must  be  measured 
and  specified.  Equipment  sensitive  to  concentrations  lower  than  1000  ppm  acetylene  must  be 
used.  Based  on  the  confused  situation  identified  in  Table  3  of  CSP,  the  reaction  order  should  be 
redetermined  and  demonstrated  to  readers.  (Towell  and  Martin^,  for  example,  claimed  second  order 
dependence,  but  their  data  indicated  an  order  much  closer  to  1.5.)  Based  on  the  literature  cited  in 
C.SP.  the  last  time  experimental  data  was  presented  showing  second  order  behavior  was  1964;  since 
then  there  have  been  many  dissimilar  results.  Considering  the  relative  sophistication  of  present-day 
experimental  techniques,  it  should  be  straightforward  to  establish  the  kinetic  order  of  acetylene 
pyrolysis. 

In  conclusion,  we  are  convinced  that  there  are  uncertainties  still  to  be  resolved  in  the  acetylene 
system.  The  trace  impurity  experiment  has  not  yet  been  performed.  An  extrapolation  of  a  rate 
expression  to  low  acetone  concentrations  suggests  a  contribution  of  10%  to  the  overall  rate  at  about 
900K.  We  feel  that  the  results  of  DACl  do  not  disprove  our  suggestion  that  even  small  concentrations 
of  impurities  can  have  a  significant  impact  on  acetylene  pyrolysis. 


3 


REFERENCES 


1.  Colket,  M  B.,  Seery,  D.J.,  and  Palmer,  H.B.,  Combust.  Flame  75:  343-366  (1989).  (CSP). 

2.  Kiefer,  J.H.  and  Von  Drasek,  W.A.,  “The  Mechanism  of  the  Homogeneous  Pyrolysis  of  Acety¬ 
lene”,  to  be  published  in  Int.  J.  Chem.  Kin.  . 

3.  Duran,  R.P.,  Amorebieta,  V.T.  and  Colussi,  A.J.,  “Radical  Sensitization  of  Acetylene  Pyroly¬ 
sis”,  to  be  published  in  Int.  J.  Chem.  Kin.  (DACl). 

4.  Colussi,  A.J.,  Combust.  Flame,  this  issue, 

5.  Duran,  R.  P.,  Amorebieta,  V.  T.  and  Colussi,  A.  J.,  J.  Am.  Chem.  Soc.  109:  3154  (1987). 
(DAC2). 

6.  Benson,  S.W,  and  O’Neal,  “Kinetic  Data  on  Gas-Phase  Unimolecular  Reactions”  ,  NSRDS-NBS 
21,  U.S.  Govt.  Printing  Office,  Wzishington,  D.C.  (1970). 

7.  Towell,  G.  D.  and  Martin,  J.  J.,  A.  /.  Ch.  E.  Journal  7:  693  (1961). 

8.  Bittner,  J.,  “A  Molecular  Beam  Mass  Spectrometer  Study  of  Fuel-Rich  and  Sooting  Benzene- 
Oxygen  Flames”,  Ph.D.  Dissertation,  Appendix  P,  Massachussetts  Institute  of  Technology. 
1981. 

9.  Ghibaudi,  E.  and  Colussi,  A.J.,  J.  Phys.  Chem.  92:  5839  (1988). 

10.  Westmoreland,  P.  R.,  “Prediction  of  Kinetics  for  C4  Species  which  form  Benzene”,  American 
Chemical  Society,  Fuel  Chemistry  Preprints  32:  480  (1987).  NOTE:  The  original  reference 
incorrectly  lists  the  heat  of  formation  of  i-C^Hs  (personal  comminication  from  P.  R.  West¬ 
moreland). 

11.  Stein,  S.  E.,  National  Bureau  of  Standards,  Gaithersburg,  MD,  Personal  Communication,  1987. 

12.  Cowperthwaite,  M.  and  Bauer,  S.  H.,  J.  Chem.  Phys.  36:  1743  (1962). 

13.  Melius,  C.,  Sandia  National  Laboratories,  Personal  Communication,  1988. 


14.  Montgomery,  J.A.,  and  Michels,  H.H.  United  Technologies  Research  Center,  Personal  Commu¬ 
nication,  1989. 


Table  I 

Reported  Thermodynamics  for  I-Buten-3-yn-2-yl  (i-C^Hs) 


AH?, 3,8 

kcal/mole 

CO 

^298 

cal/mole/K 

300 

500 

—  cal/mole/K 
800  1000 

1500 

cm® /mole/sec 

Ref. 

110.8 

63.7 

16.3 

21.9 

27.8 

29.8 

33.1 

0.37 

8 

115.2 

66.1 

16.9 

22.1 

27.0 

29.3 

33.0 

0.11 

9 

114.5 

68.3 

17.5 

22.5 

27.4 

29.6 

33.2 

0.56 

10 

116.1 

68.1 

16.8 

22.2 

27.1 

29.3 

32.8 

0.19 

11 

102.0 

65.4 

17.4 

22.4 

27.1 

29.2 

32.5 

126 

12 

110.0 

71.6 

17.4 

22.4 

27.1 

29.2 

32.5 

33.8 

13 

*  assuming  a  k_ia  =  5xl0'^  cm®/mole  sec 

Note:  The  low  heat  of  formation  calculated  by  Ref.  1 1  is  probably  in  error,  but  is  the  value  included 
in  many  thermodynamic  compilations.  The  relatively  high  value  of  entropy  determined  by 
Melius'^  has  recently  been  supported  by  additional  calculations*^. 
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Tie  ‘’.r.'tic?  cf  r;Trlytis  of  pyridire  dilo.ce  1  r-  areor  have  ''eep 
■;  r  =  ~  ;  "2  I  e-oc  i  se  shoe'"'  tiif'e,  csirr  capi’larv  ceforr  CC- 
tO't'ether  wirh  rC/'-'S  ard  FTIF  spectroscopy  for  product  de  terri  r.a  t  i  or  , 
over  t^'e  terperature  ranee  !3d'0-lF00f'  and  total  pressures  of  F-ll 
ato.  Ac  cHe  lower  end  cf  the  studied  cerperatures .  cyanoacetvlene  was 
found  to  be  the  principal  nitropen-containine  product.  At  elevated 
cerperatures  hydrogen  cyanide  predorina ted .  Other  ra jor  products  were 
acetylene  and  hydrogen. 

Therr.ocherical  escirates  of  the  isoreric  cyclic  pyridyls  produced 
in  the  pyrolysis  indicate  that  the  ortho-isomer  is  uniaue  in  being  able 
to  undergo  facile  cleavage  to  an  open-chain  cyano-rad ica 1  from  which 
cyanoacetylene  is  produced.  Several  sources  of  HC.V  were  identified  in 
the  system.  The  r-  and  p-pyridyls  ray  eliminate  HCN  in  a  molecular 
process-  An  important  source  of  HCN  at  high  temperatures  is  the 
addition  of  H  atoms  to  cyano-compounds .  especially  cyanoacetylene,  but 
also  acetonitrile  and  acrylonitrile  which  are  produced  in  the  pyrolysis. 

The  pyrolysis  is  a  chain  process  initiated  principally  by  C-H  bond 
fission  tc  form  o-pyridyl-  A  5f-scep  reaction  model  is  nresented  and 
shown  to  substantially  fit  the  observed  profiles  of  the  major  product 
species.  From  this  model  we  derive  a  value  for  the  rate  constant  of  the 
principal  initiation  reaction 

-  o-C^H^.v  +  H  (1) 

Qf  *  j^qIS. 9*0.1  jjxp  (-98*3^031  ~o  1  ~  'PT )  s“  ‘  between  1300-1800K  and  at 
a  total  pressure  of  about  1C  at-. 


I  r,  c  r  0  u  c  t  i  o  r. 


i'caw  ruels  cue';  as  coal  and  coal-der i  vpd  liquids  contain 
cor. '^ic'era^'Ic-  arouncs  of  cberically  bound  nitropen  and  the  corbustion  of 
these  fuels  can.  lead  to  significant  arounts  of  fuel-derived 
''uch  of  the  nitrogen  originally  bound  in  the  fuel  is  in  the  fort?  of 
aroratic  hete rocyc 1 i c  structures  containing  pyridine  and  pyrrole  ring 
systers^'^.  Lnder  pyrolytic  conditions  these  heterocycles  tray  fort? 
nitrogen  precursors  of  NO^  so  that  the  rates  of  pyrolysis  and  product 
forration  frot?  representative  nitrogen  heterocycles  will  determine  NOj^ 
formation  in  the  combustion  of  heavy  fuels. 

— P 

The  thermal  decom.pos i t ion  of  pyridine  has  been  studied  previously 
in  flow  systems  between  1100-1400K.  Acetylene  and  hydrogen  cyanide  have 
been  found  to  be  major  products.  Other  nitrogen-containing  products  have 
been  identified  as  acrylonitrile,  benzoni tr i le ,  quinoline  and 
acetonitrile,  although  dispute  exists  over  the  observation  of 
■  ^tonitrile. 

Axworthy  et  al.^  found  that  the  decomposition  of  pyridine  followed 
first  order  kinetics  with  an  overall  rate  constant  of  3.8  x  10^^  exp  (-70 
kcal  mol^VRIjs"^  between  1100-1200K.  Houser  et  al.^  found  that  the 
decomposition  kinetics  were  complex  and  did  not  follow  first  order 
kinetics  exactly.  However,  they  derived  an  approximate  first  order  rate 
constant  of  1  x  10^^  exp  (-75  kcal  rol~VRT)s~^,  about  an  order  of 
magnitude  lower  than  that  of  Axworthy  et  al. 

Recently,  Kern  et  al^  carried  cut  a  shock  tube  study  of  the  pyrolysis 
of  pyridine  over  the  temperature  var^c-  of  I7CO-24CO  K  and  at  pressures 


e  t;  We  e  n  _  . 


3  cCr  usirp  hoc'f  laser  scblieren  dens  i  ton-e  t  ry  and  Cire-of- 


rass  s  pec  C  ror  e  t  r  y  as  diaproscic  tods.  Tl'ey  proposed  a  chain 
■'echanisr  wic'r  ’irsc-crder  initiation  involving  a  rinc  hydroeen  fission 

C^H.N  -  C,H,N'  H  (la) 

5  0  0  i 

wit!'  a  specific  pressure-dependent  rate  constant  of  1.08x10^2  exp(-74.5 
!<cal  rol  -/RT)  s  In  their  analysis  the  represented  a  weighted 

average  of  the  isoreric  pyridyl  radicals.  Froir  a  RPKy  extrapolation  of 
their  rate  data,  Kern  et  al  derived  an  average  value  of  76±10kcal  rol~^ 
for  the  heat  of  forration  of  the  pyridyl  radicals,  a  value  higher  than 
that  i-plied  by  the  earlier  flow  studies^*®’ 

In  a  recent  shock  tube  study  Leidreiter  and  Wagner^®  iteasured  the 
rate  of  decomposition  of  pyridine  at  temperatures  between  1700-2000K  and 
at  total  densities  of  2  x  10  -3  x  10  mol  cm  ,  using  time  resolved 

absorption  measurements  at  279  nm.  They  did  not  determine  products  but 

from  fall-off  curves  they  were  able  to  extrapolate  their  rate  data  to 
obtain  the  limiting  high  pressure  rate  constant  for  disappearance  of 
pyridine  of  k«  =  10^^-2  exp  (-100  kcal  mol'VRT)  s“^. 

The  lower  reported  values  of  the  activation  energy  of  the  initiation 

step  (/0-75  kcal  mol“^)  imply  that  the  C-H  bond  energy  In  pyridine  is 
considerably  lower  than  that  in  benzene  (assuming  that  the  initiation  step 
is  C-H  fission).  Also,  Kern  et  al.'^  detected  only  one  nitrogen-containing 
product  -  hydrogen  cyanide  -  although  in  flow  studies  several  other 
nitrogen  products  were  observed.  Thus,  the  present  shock  tube  kinetic 
study  of  pyridine  pyrolysis,  using  capillary  column  gas  chromatography  and 
chromatography/mass  spectrometry  for  detection,  is  oriented  towards 


f^ecsilec  proci.'ct  i  den  c  i  f  ica  c  iop.  and  elucidation  of  the  pyrolysis  ir-echanisr 

i'xper  irenta  1 

Pyrciyses  were  carried  out  on  two  sinple-pulse  shock  tubes,  one  of 
diareter  3.8  c~  at  L'TRC  which  was  pritrarily  used  for  kinetic  studies,  the 
other,  at  the  Lniversity  of  Sydney  (diameter  7.6  cm),  was  used  principally 
for  product  identification  in  conjunction  with  FTIP  and  a  CC/VS  at 
CSIRG.  Eoth  shock  tubes  are  basically  similar;  the  I'TRC  tube  and  on-line 
capillary  CC.  analytical  system  have  been  described  previous  ly^  ^ . 

Pyridine  samples  were  purchased  from  Aldrich  (HPLC  grade)  and  from 
Merck  (A.P.  grade).  Further  purification  by  several  bulh-to-hulb 
distillations  produced  pyridine  samples  of  purity  exceeding  99. 9?^  (by  CC 
with  CH^CN  as  the  only  detectable  impurity).  No  difference  between  the 
two  sources  of  pyridine  was  observed. 

For  the  kinetic  studies,  mixtures  of  0.7  and  0.15%  pyridine  in  argon 
were  prepared  and  stored  in  glass  bulbs.  Pressures  and  t.emperacures 
behind  the  reflected  shock  were  computed  from  the  measured  incident  shock 
velocity.  Residence  times  and  quenching  rates  by  the  rarefaction  wave 
were  measured  from  pressure  profiles  recorded  by  Kistler  gauges. 

Reflected  shocked  gas  temperatures  ranged  from  1200-1800K  and  total 
pressures  from  7-11  atm.  Residence  times  at  uniform  temperatures  behind 
the  reflected  shock  front  ranged  from  450-600  us.  For  the  product 
analyses,  three  separate  capillary  CC  columns  were  ured.  The  L’TRC  shock 
tube  was  eouipped  with  a  Chrompack  CP  Sil  5  capillary  column.  The  Sydney 
shock  tube  employed  a  SCE  BP-1  (5um  bonded  phase  d imethysi loxane)  wide- 


bore  colurn  vhereas  che  CSIRO  GC/>'S  was  eauipped  with  an  SGE  BP-1  0.5  un 
bonded  phase. 

'.'lietc  corrercial  samples  were  available,  product  identification  ana 
calibration  were  rade  from  capillary  CC  retention  times  and  areas  by 
direct  corpari^or.  with  standards-  HCK  was  purchased  as  a  standard  mixture 
in  N'2  fi'o'"  Scott  Specialty  Cases.  Acetonitrile,  acrylonitrile  and 
benzonitrile  (all  >  997  purity)  were  obtained  from  Aldrich. 
hydrocarbons  were  calibrated  gas  standards  prepared  and  analyzed  by  Scott 
Environmental  Laboratories. 

Cyanoacety lene  (HCCCN)  was  identified  in  the  product  gases  by  FTIR 
absorption  spectroscopy  in  a  7m  cell  (via  its  strong  CCH  bend  at  663. A  cm 
and  the  C.V  stretch  at  2272  cm~^).  Product  samples  from  the  Sydney  shock 
tube  were  collected  in  gas  burettes  and  transported  to  CSIRO  for  GC/MS 
analysis.  The  CC  peak  assigned  to  cyanoacetylene  was  confirmed  not  only 
by  the  observation  of  the  parent  peak  at  51  a.m.u.,  but  also  through  the 
observation  of  its  fragmentation  pattern  whose  intensities  were  in 

1 2 

excellent  agreement  with  those  reported  by  Dibeler  et  al 

Other  products  confirmed  both  by  FTIR  and  CC/MS  were  hydrogen 
cyanide,  vinylacetylene  and  diacetylene.  In  addition,  triacetylene  (C£,H2) 
and  cyanodiacetylene  were  positively  identified  by  GC/MS. 

Both  hydrocarbons  and  nitrogen-containing  compounds  were  detected  by 
FID;  the  sensitivity  of  detection  of  HCN,  however,  was  only  about  one 
tenth  that  of  C.,H2.  permitting  measurements  of  HCN  to  be  made  with 
adequate  precision  above  about  30  ppm.  Hydrogen  was  detected  by  TCD. 


Tl'e  capillriry  c c  1  iirr s arp  1  i ne  svster  exhibitec'  a  "rerory  effect"  for 
rV-e  r  :  C  roper.-ccn  ta  i  p  i  ne  corpounds,  especially  pvridine,  ben  zon  i  t  r  i  1  e  , 

•ice  ton  i  c  r  i  le  and  acrylonitrile.  Repeatability  of  injections  of  these 
ccnpounds  was  only  ±  Pyridine  mixtures  in  argon  slowly  decreased  in 

concentration  with  storage,  requiring  frequent  recalibrations. 

Adsorption  and/or  polymerization  of  the  product  cyano-cotrpounds  on  to 
the  walls  of  the  gas  burettes  was  a  problem  leading  to  significant  losses 
in  sensitivity  of  these  products  in  transportation  of  samples  for  GC/MS 
analyses  . 

Results 

Yields  of  all  products  of  significance  (scaled  by  the  initial 
concentration  of  pyridine)  are  shown  in  Figure  1  as  a  function  of 
temperature  for  the  0.77-  mixture  of  pyridine  in  argon.  The  variation  of 
pyridine  concentration  with  temperature  is  also  given  in  Figure  1. 

The  principal  hydrocarbon  product  at  all  temperatures  and  for  both 
series  of  mixtures  was  acetylene.  Hydrogen  was  an  important  product  at 
all  tem.peratures .  The  C2H2/H2  ratio  varied  from  about  1/3  at  1300K  to 
about  1.5/1  at  1600K.  At  the  lower  end  of  the  studied  temperature  range 
cyanoacetylene  was  the  principal  nitrogen-containing  product  being  higher 
in  yield  than  HCN  by  about  a  factor  of  2.  From  FTIR  difference  spectra  of 
the  products  from  successive  runs  in  the  lower  range  of  temperatures  it 
was  shown  that  the  rates  of  increase  with  temperature  of  both  acetylene 
and  cyanoacetylene  were  approximately  equal.  With  increasing  temperature 
HCN  yields  overtook  those  of  HCCCN  such  that  the  yield  of  HCN  is  about  ^ 
tires  that  of  HCCCN  at  IPOOK.  Nuch  smaller  yields  of  acetonitrile  and 


a  c  r  y  1  on  i  t:  r  i  It  were  recorc'et^  at  all  teroe  ra  t;  ur  e  ?  .  An  additional  nroduct 
-  i  s  '^erected  and  vas  assured  to  he  a  r;-'ano-cor  ta  i  n  i  rg  species  due  both  to 
: e  characteristic  tailing  of  the  GC  peak  t^at  this  class  of  compounds 
exhibits  on  the  capillary  column  and  also  to  Che  lack  of  appearance  of 
tlis  peak  during  pvrolysis  of  pure  hydrocarbons.  This  product  peaked  in 
yield  at  about  15CCK  and  was  present  even  at  the  lowest  temperatures  at 
which  decoitposi c ion  could  be  detected.  From  its  retention  time,  eluting 

just  after  C^H2,  which,  in  turn,  eluted  soon  after  HC5‘i  it  been 

tentatively  assigned  as  cyancvinylacetylene  ( S-e thynylacryloni tri le : - 
(HC^-CK=CHCN) .  Actem.pts  to  confirm  this  assignment  by  GC/MS  were 
unsuccessful,  however,  as  the  level  of  this  peak  after  transportation  of 
the  product  gases  for  analysis  had  dropped  below  sensitivity  limits. 
Cyancvinylacetylene  is  known  to  polymerize  rapidly^^  above  -30°C. 

Other  products  included  vinylacetylene  (C4F4)  -  important  at  the 
lower  temperatures,  diacetylene  and  triacetylene  -  important  over  the 
entire  temperature  range,  methane,  ethylene  and  benzene.  There  were  also 
significiant  yields  of  benzoni tr i le ,  allene  and/or  propyne  (unresolved  by 
the  capillary  columns)  together  with  traces  of  propene,  butadiene,  styrene 
and  phenylacetylene .  Cyanod iacetylene  (HC^N’)  yields  only  became 
significant  at  the  higher  temperatures  of  this  study. 

When  results  from  the  two  series  of  initial  pyridine  concentrations 
(0.7  and  0.157()  were  compared,  it  was  found  that  the  overall  decomposition 
of  oyridine  and  formation  of  major  products  HCN,  C2H2,  HCCCN  and  H2  were 
essentially  independent  of  initial  pyridine  concentration.  This 
independence  upon  initial  pyridine  concentration  was  also  observed  for 
products  C^H^,  CH4  and  acrylonitrile.  Thus,  overall  first  order  kinetics 


are  observer  fcr  pyridine  dccor-ros  i  t  ion. ,  and  for  tbe  for~acion  of  cbe 
ra’^r  as  veil  as  •'anv  of  tbe  rinor  species.  Tbe  first  order  be^^avior  is 
consistent  vitb  that  previously  observed-'  Because  of  the  first  order 
bei-avicur  it  is  t''erefore  convenient  in  reporting  species  data  to  present 
nor-alized  concent ra t i cns  of  reactant  and  products  in  Figures  1  and  3-b. 

DISCL'SSION 

Tberr.ocherristry  of  the  Pyridine  System 


The  reaction  n-echanisir  developed  to  interpret  the  above  results 
involves  a  chain  reaction  initiated  by  C-H  bond  scission,  as  suggested 
previously.  Subseauent  reactions  are  dependent  on  the  thertnochemist ry  of 
the  cyclic  pyridyl  radicals  produced  by  the  scission.  The  enthalpies  of 
forration  of  the  pyridyls  were  estirated  under  the  assumption  that  the  C-H 
bond  energy  in  pyridine  was  equal  to  that  in  benzene  together  with  a 
stabilization  energy^^  of  1.5  kcal  mol”^  for  the  o-  and  p-pyridyls.  Using 
the  assumption  of  equality  of  C-H  bond  energies  of  pyridine  and  benzene, 
Dewar  and  co-worker s have  calculated  enthalpies  of  atomization  of 
pyridine  and  other  molecules  containing  the  pyridine  ring,  in  very  good 
agreement  with  experiment  by  SCF  MO  techniques.  Heat  capacities  and 
entropies  of  the  pyridyls  have  been  estimated  by  group  additivity 
methods^^  by  assuming  that  the  contribution  of  the  Cg  group  to  the  pyridyl 
thermochemistry  was  equal  to  that  in  phenyl. 

All  pyridyl  radicals  may  undergo  thermal  ring  fissure.  However,  the 
o-pyridyl  is  unique  in  its  ability  to  produce  directly  an  o''en-chain 
cyano-rad ical  (I).  Figure  2  shows  all  possible  open-chain  "linear" 
radicals  l-C^H/.v  arising  from  simple  ring  cleavage.  As  the  values  in 
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Table  I  indicate,  the  cyano-rad ica 1  (I)  should  be  overwhelmingly  more 
stable  that  any  ocher  open-chain  radical  produced  by  ring  fissure  of 
pyr idyls . 

Because  of  the  lack  of  thermochemica  1  data  for  the  pyridyls  we  have 

had  to  resort  to  approximate  group  additivity  methods  to  arrive  at  the 

estimates  of  the  thermochemical  parameters  given  in  Table  I.  To  estimate 

thermochemical  parameters  for  radical  (I),  thermochemical  data  for  group 

C^-(H)(CN)  were  obtained  by  additivity  methods^'^  from  the  acrylonitrile 

thermochemistry  together  with  data  for  n-C^H^  (Ref.  16).  Estimates  for 

radical  (II)  were  calculated  from  the  thermochemistry  of  CH=NH,  l-CgH5  and 

C2H2  (Ref.  16).  Thermochemical  data  for  CH=NH  were  calculated  by  Melius 

using  the  BAC-MP4  technique^^.  The  enthalpy  of  formation  of  l-Cg^Hj  was 
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estimated  by  additivity  methods  and  the  value  of  A  »  139  kcal 

mol~^  used  in  the  present  work  compares  well  with  a  recent  determination 
by  Braun-Unkhof f  et  al^^  of  144±3kcal  mol~^.  group  data  (required 

for  radical  (III))  were  estimated  by  comparing  group  differences  between 
_^C=CH  anc^C=CH2  witl^C=N  and  ^C=NH  where  Nj-(H)  data  are  available  from 
CH2=NH.  The  enthalpy  of  formation  of  CH2=NH  has  been  measured 
experimentally^*^.  Heat  capacities  and  entropies  are  obtained  from 
spectroscopy^^. 

Thermochemical  data  for  HCCCN,  and  H2C=CHCN  were  computed  by 

statistical  methods  using  molecular  constants  from  Refs.  22-24.  Other 
relevant  group  data  are  available  in  the  thermochemistry  literature . 


^eac  c_^n _ >!e_(^iar  i  sr 

Easpd  on  the  assurpc  ion  of  tl'e  stahilicy  of  the  open-chain  cvano- 
radical  (i),  the  dorrinant  rrechanistr’  in  the  pyrolysis  of  pyridine  in  the 
lower  end  of  the  studied  terperature  range  of  1300-1S00K  is  postulated  to 

be  (Reaction  nurberinp  as  in  Table  II): 

°  -1 
^H3oo/l<cal  mol 


C^H^N  -  o-C^H^N  +  H  109  ±  3  (1) 

o  -  -  HC=CH-CH=CHCN  (I)  22  ±  6  (3) 

HC=CH-CH=CHCN  ♦  +  HC=CHCN  37  ±  6  (4) 

HC  =  CHCi;  (+M)  -  HCCCN  +  H  C+M)  40  ±  3  (5) 

H  +  C^H^N  -  0-.  P-,  r-CjH^N  +  (6,7) 


This  is  the  major  low  temperature  chain  mechanism  although  a  minor 
chain  process  may  arise  via  fission  of  a  secondary  H  from  (I),  producing 
cyanovinylacetylene  (C\'A). 

HC=CH-CH=CHCN  -  HCHC-CH=CHCN  +  H  =  38  ±  3  kcal  mol"^  (8) 

The  HCCHCN  radical  was  found  to  decompose  principally  via  Reaction 
(5).  Other  possibilities  Include 

HCCHCN  *  +  CN  =  61  Veal  mol"^ 


KCCHCN  -  C\H  +  HCN  ^  ^ 

but  the  larger  endothermic ities  of  these  reactions  make  them  less  likely 
at  least  at  the  lower  temperatures  of  this  study. 

The  above  short  chain  mechanism  is  in  accord  with  the  experimental 
observations  that  HCCCN'  is  the  principal  nitrogen-containing  product  at 
low  temperatures,  and  of  the  equivalence  in  production  rates  of  C2H2  and 
HCCCN.  It  also  gives  an  explanation  for  the  product  peak  (assigned  as 
CVA)  which  appears  at  low  temperatures. 

The  thermochemistry  of  the  m-  and  p-pyridyls  suggests  that  these 
radicals  will  be  stable  with  respect  to  simple  ring  cleavage,  hence  the 
reverse  of  Reaction  (2) 

p-,m-C.H,N  +  H  ♦  C,.H,,N  (-2) 

5  «  5  5 

terminates  radicals.  However,  both  of  these  isomers  may  also  undergo 
elimination  of  HCN  to  produce  the  n-C^H3  radical  (Reaction 
(9)).  There  are  other  possible  sources  of  HCN.  Open-chain  radical  (I) 
could  eliminate  HCN  in  a  molecular  process  involving  a  cyclic  transition 
state  (Reaction  (10)).  Although  we  are  unaware  of  any  previous  studies  of 
HCN  elimination  from  cyclic  radicals,  Braun-Unkhof f  et  al^^  have  analysed 
the  decomposition  of  benzyl  radicals  by  invoking  a  rearrangement  to  a 
bicyclic  O^H^  radical  which  can  molecularly  eliminate  C2H2  in  an 
analogous  process  to  our  HCN  elimination. 

At  elevated  temperatures  where  high  concentrations  of  H  atoms  exists, 
an  important  source  of  HCN  arises  from  the  addition  of  H  atoms  to  cyano- 


compounds.  The  most  important  of  these  is  reaction  (-27)  which  leads  to 
the  decay  of  HCCCN  at  high  temperatures.  The  abstraction  reaction  of  CN 
with  H2  also  plays  a  significant  role  at  elevated  temperatures. 

The  detailed  reaction  mechanism  postulated  to  simulate  the  yields  of 
the  major  products  is  presented  in  Table  II.  In  addition  to  the  nitrogen 
chemistry,  the  mechanism  includes  acetylene  and  benzene  formation 
reactions  as  subsets.  For  many  of  the  reactions  involving  nitrogen  - 
containing  species,  few  or  no  rate  constant  evaluations  are  available  in 
the  literature.  For  these  reactions  estimates  were  made  using  rate 
constants  from  analogous  reactions  or  thermochemistry.  Based  on  reaction 
path  analyses  and  sensitivity  analyses,  selected  rate  constants  were 
adjusted  to  fit  the  experimental  data.  These  latter  rate  constants  are 
denoted  "PW"  in  Table  II. 

For  simplification,  only  one  radical  is  considered  and  this  is 

0 

taken  to  be  the  n-  radical  whose  is  assumed  to  be  125  kcal 

mol”^.  Allene/propyne  is  c  ’ 'idered  to  be  propyne,  the  more 
thermodynamically  favoured  form  of  at  high  temperatures.  The 

reaction  model  includes,  in  addition  to  the  short  chain  mechanism  referred 
to  previously,  secondary  reactions  of  the  product  species  and  radical- 
radical  reactions.  There  appear  to  be  strong  parallels  between  the 
addition  reactions  of  acetylene  and  those  of  cyanoacetylene .  Thus  rate 
constants  for  radical  additions  to  HCCCN  have  been  chosen  by  analogy  with 
literature  constants  for  additions  to  C2H2.  Also,  acrylonitrile  should  be 
able  to  undergo  2,2  and  3,3  elim.lnations  of  HCN  and  H2  respectively, 
analogous  to  the  1,1  elimination  of  H2  and  vinylldene  formation  reactions 
of  ethylene  and  vlnylacetylene^^.  Integration  of  the  kinetic  equations 


was  made  using  the  CHEMKIN  package-7a^  LSODE27b  g^d  the  SANDIA  shock  tube 


code* 


which  was  rodirief*  to  account  for  the  effects  of  Quenching.  Rate 


Ltf  production  and  sensitivity  analyses  ^'ere  carried  out  by  the  FF.SKiy 

^  -B 

code 


Corparison  of  the  kinetic  rodel  with  experipent  is  given  in  Figures 
3-6  for  the  rajor  products  (i-e-,  species  whose  yields  are  >0.1  over  the 
studied  temperature  range).  The  model  adequately  describes  the 
temperature  dependence  of  the  major  products  and  minor  products  also  from 
both  series  of  runs.  The  model  also  reproduces  the  experimental  pyridine 
profiles  up  to  about  80%  decomposition.  Although  about  5%  of  the  initial 
pyridine  remains  unreacted  due  to  boundary  layer  effects  at  1800K,  this 
has  little  effect  on  the  profiles  of  the  major  species.  For  example,  a 
decrease  of  about  5%  in  the  HCN'  concentration  at  1800K  is  well  within  our 
experimental  error. 

It  should  be  noted  that  the  mechanism  of  Table  II  contains  certain 
reactions  postulated  to  simulate  minor  products'  formation.  Rate 
constants  were  determined  based  on  best  fits  to  experimental  data  and/or 
estimates  as  necessary.  In  particular,  CF4.  allene/propyne  and 
acetonitrile  were  observed  experimentally  in  appreciable  yields  at  quite 
low  temperatures,  implying  their  formation  through  reactions  of  low 
activation  energy.  To  simulate  the  formation  of  these  minor  products,  a 
route  of  low  activation  energy  via  H  addition  to  pyridine  has  been 
postulated.  Reactions  (30)  and  (31)  are  actually  a  simplification  for  a 
com.plex  addition  and  ring  contraction  mechanism:- 


•  ^ 


The  .Tiechyl pyr rolyl  radical  (IV)  produced  in  the  above  scheme,  should 
possess  considerable  resonance  stability  as  it  is  analogous  to  the 
re  thy Icyc lopent ad ienyl  radical-  The  group  additivity  estimate  for 
(IV)  “  46  kcal  mol~^  (Table  I)-  The  methylpyrrolyl  may  thus  be  formed  in 
a  chemically  activated  state-  It  is  further  suggested  that  (IV)  may 
undergo  rapid  unimolecular  reaction  to  CH3CN  and  propargyl  radical 
(C3H3).  Methane  can  then  arise  from  methyl  radicals  produced  by  reaction 
of  CH3CN  with  H  (Reaction  29).  Propargyl  radicals  are  an  Important  source 

of  the  C3H4.  The  proposed  reaction  scheme  to  (IV)  is  analogous  to  an 
addition  and  ring  contraction  mechanism  postulated  by  Ritter,  Bozzelll  and 
Dean^^  for  H  addition  to  benzene  andmight  be  regarded  as  controversial. 
Nevertheless,  these  and  other  reactions  of  the  mechanism  postulated  to 
simulate  formation  of  the  minor  species  have  only  a  slight  effect  on  the 
major  products'  profiles,  and  have  been  included  solely  for  completeness. 

A  brief  selection  of  results  of  sensitivity  analyses  of  the  reaction 
mechanism  is  given  in  Table  III.  This  Table  contains  normalized 
sensitivity  coefficients  evaluated  at  1600K  and  at  100  us  for  the  mixture 
of  0.7Z  pyridine  in  argon.  The  accuracy  and  usefulness  of  a  sensitivity 
analysis  is  dependent  strongly  on  the  validity  of  the  selected  mechanism 
as  well  as  the  specific  rate  constants.  Considering  that  this  detailed 
chemical  kinetic  model  of  pyridine  decomposition  is  the  first  to  describe 
the  formation  of  cyanoaceftjiene  and  several  other  important  nitrogen- 
containing  products  and  intermediates,  the  results  in  Table  III  should  be 
considered  qualitative  rather  than  quantitative.  This  fact  is 
particularly  true  for  portions  of  the  model  that  are  speculative.  For  the 
main  reaction  sequence  which  is  considered  less  speculative,  the 
sensitivity  analysis  indicates  that  the  reactions  most  important  to  the 
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decay  of  the  reactant  and  formation  of  major  products  are  Reactions 
1,4, 5,6  and  7.  The  low  sensitivity  to  ^3  (the  sensitivity  coefficient  for 
all  the  species  listed  in  Table  III  to  this  rate  constant  was  less  than 
10“^)  indicates  that,  at  least  for  the  value  of  k3  estimated  in  this  work. 
Reaction  3  is  not  a  rate  limiting  step  to  thermal  decomposition  of 
pyridine;  rather,  it  is  the  production  of  the  o-pyridyl  radical  via 
Reaction  6  (at  1600K)  which  is  rate  limiting.  At  lower  temperatures 
(1300-1400K) ,  the  value  of  the  sensitivity  coefficient  for  the  pyridine 
concentration  with  respect  to  kj^  (at  100  microseconds)  is  a  factor  of  four 
greater  than  for  that  with  respect  to  any  other  rate  constant.  Under 
these  conditions,  Reaction  1  is  a  rate-limiting  step  for  the  formation  of 
o-pyrldyl  as  well  as  for  the  overall  decomposition  of  pyridine. 

Our  optimized  value  of  the  rate  constant,  ■  10^^’^  ^  exp 
(-98±3kcal  mol”VRT)&“^,  is  in  excellent  agreement  at  1700K  with  the 
limiting  high  pressure  rate  constant  obtained  by  Leldreiter  and  Wagner^® 
for  disappearance  of  pyridine.  Our  value  differs  from  theirs  by  less  than 
5X  at  this  temperature,  and  both  sets  of  Arrhenius  parameters  agree  well 
within  experimental  error.  Thus,  based  on  Leldreiter  and  Wagner’s^®  fall- 
off  calculations  we  may  assume  that  under  our  experimental  conditions  (*10 
atm  pressure)  the  unlmolecular  decomposition  of  pyridine  was  caking  place 
at  or  very  close  to  the  high  pressure  limit  and  that  Leldreiter  and 
Wagner's^®  rate  constant  for  disappearance  of  pyridine  may  be  equated  with 
our  value  of  k^,  rate  constant  for  the  fission  of  the  ortho-hydrogen  of 
pyridine. 

The  reaction  mechanism  of  Table  II  contains  several  other  pressure- 
depeodei.c  rate  constar..s.  Some  of  these  reactions  may  be  in  the  fail-off 
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region  and  Che  rate  constants  for  these  reactions  may  therefore  be 
specific  for  the  present  conditions. 

Our  values  of  and  k7  (Table  II)  imply  experimental  heats  of 
formation  for  o-pyridyl  radical  of  ®  82±6  kcal  mol  ^  and  for  m-  and 

p-pyridyls  of  86  ±  6  kcal  mol“^.  The  latter  value  just  encompasses  the 
chermochemical  estimate  (Table  I)  within  its  error  limits,  however,  the 
difference  in  heats  of  formation  of  o-  and  the  m-,  p-pyridyls  is  larger 
than  the  thermochemical  estimate^®.  An  ortho-stabilization  energy  as 
large  as  4  kcal  is,  however,  predicted  by  ab  initio  molecular  orbital 
calculations  on  o-pyridyl^*^  in  which  there  is  appreciable  stabilizing 
interaction  between  the  nitrogen  non-bonding  orbital  and  the  singly 
occupied  molecular  orbital  on  the  ortho-carbon.  This  interaction  is  not 
possible  for  the  m-  and  p-radicals. 

Our  values  of  for  the  pyridyls  are  also  larger  but  within 

uncertainties  of  76±lO  kcal  mol”^  obtained  as  an  average  value  for  the 
pyridyls  by  Kern  et  al^  using  an  RRKM  extrapolation  of  their  rate  data. 
However,  in  view  of  the  possible  effects  of  trace  Impurities, 
uncertainties  in  single  pulse  shock  tube  measurements  at  large  extent  of 
conversion  and  the  additional  possibility  that  the  reaction  might  still  be 
in  the  fall-off  region,  we  cannot  yet  positively  reject  the  thermochemical 
estimates  of  200  pyridyls  in  favor  of  our  experimental  values. 

Our  kinetic  model  when  applied  directly  to  simulate  the  TOF  profiles 
of  Kern  et  al.^  overpredicts  the  rates  of  decomposition  of  pyridine  and 
formation  of  C7H2  although  its  prediction  of  HCN  rates  is  in  agreement. 
Since  our  model  was  developed  for  a  higher  range  of  pressures  and  in  view 
of  the  likelihood  that  the  unlmolecular  initiation  reactions  (1)  and  (2) 
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are  likely  to  be  well  into  the  fall-off  region  at  their  pressures,  the 

Q 

disagreement  is  not  surprising.  The  failure  of  Kern,  et  al  to  observe 
cyanoacetylene  amongst  their  products  is  probably  not  inconsistent  with 
our  results  -  we  estimate  that  at  their  highest  temperature  (2092K)  HCCCN 
should  be  only  one  fifth  of  the  concentration  of  HCN.  At  the  lowest 
temperature  of  their  studies  (1703K)  HCCCN  and  HCN  should  be  of  comparable 
concentration.  However,  sensitivity  problems  might  then  obscure  mass 
spectral  detection  of  HCCCN,  especially  as  pyridine  itself  has  a  strong 
fragment  peak  at  51  A.M.U.  (Ref.  31). 

Sooting  Tendency  of  Pyridine 

In  many  previous  studies  (see  references  in  Leidreiter  and  Wagner^*^) 
the  sooting  tendencies  of  pyridine  in  flames  or  shock  tubes  have  been 
found  to  be  substantially  less  than  that  of  other  aromatics  (e.g. 
benzene).  The  results  of  this  investigation  are  consistent  with  the 
previous  studies  since  (except  for  benzonlti ile)  very  few. species  with 
molecular  weights  higher  than  pyridine  were  observed.  Negligible 
quinoline  and  naphthalene  were  detected.  Furthermore,  the  carbon  balance 
was  reasonable  (i;lSZ)  considering  the  difficulty  in  measuring  accurately 
the  concentrations  of  pyridine  and  other  nitrogenous  species.  This 
situation  differs  dramatically  from  the  pyrolyses  of  toluene-^^  and 
benzene^®  which  exhibited  carbon  losses  of  70Z  or  more  (pre-umsbly  to 
polyaromatic  species  and/or  soot).  In  a'^dition  several  multi-ringed 
species  were  ob^rved,  e.g.,  indene  and  naphthalene. 

Explanations  for  the  different  sooting  characteristic  of  pyridine 
usually  include  the  fact  that  with  an  N-atom  In  the  ring  large 
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per icondensed  polyaromatic  compounds  cannot  be  formed.  Based  on  the 
present  experimental  results,  it  appears  that  it  is  also  difficult 
(although  not  impossible)  to  grow  a  second  ring.  We  offer  two  possible 
explanations  for  this  phenomenon.  The  most  likely  explanation  is  the 
short  lifetime  of  the  pyridyl  radical  which  is  less  than  a  nanosecond  at 
temperatures  of  this  study.  Our  rate  constant  used  in  this  half-life 
calculation,  is  based  on  the  endothermlcity  of  Reaction  3  and  a  selection 
of  10^^  for  the  A-factor.  Although  the  model  calculations  were 
insensitive  to  this  rate  constant,  we  estimate  a  potential  uncertainty  of 
up  to  one  and  one-half  orders  of  magnitude.  In  comparison,  the  half-life 
of  the  phenyl  radical  is  about  50  microseconds  (or  500  microseconds  using 
the  recent  rate  data  of  Braun-Unkjof f  et  al^^  at  1500K.  Phenyl  radical, 
in  turn,  is  believed  to  be  a  critical  intermediate  leading  to  the  growth 
of  polyaromatic  hydrocarbons.  The  calculated  life-time  of  o-pyrldyl  is 
unlikely  to  be  in  error  by  as  much  as  four  orders  of  magnitude  and  so  we 
conclude  that  its  relatively  short  lifetime  prevents  it  from  playing  a 
critical  role  in  the  growth  of  two-ringed  species.  (The  other  pyridyls 
will  attain  larger  concentrations  but  due  to  rapid  Interchange  with  o- 
pyridyl,  their  concentrations  will  also  be  low.)  A  second  more 
speculative  explanation  relates  to  the  relative  reactivity  of  the  carbon 
atoms  on  pyridine.  We  speculate  that  the  carbon  atom  with  the  lowest 
reactivity  is  located  at  the  raeta  position.  Analogous  to  the  growth 
mechanism  for  polyaromatic  hydrocarbons,  acetylene  will  add  to  pyridyl  but 
specifically  to  o-pyridyl  or  p-pyridyl.  An  addition  of  a  second  acetylene 
and  subsequent  ring  closure  will  be  inhibited  since  external  bonding  to 
the  N-atom  will  result  in  the  loss  of  aromaticity  and  since  the  meta 
carbon  is  relatively  unreactive.  Growth  to  form  a  two-ringed  compound 
will  therefore  be  suppressed  relative  to  that  for  a  pure  hydrocarbon. 
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TABLE  I:  Thermochemical  Parameters  for  the  Pyridine  System 


0 

Species  ^^f,300  ^ 

(kcal  mol  (cal  mol”*'  K”*)  (cal  mol”*  K  *) 


300 

500 

1000 

1500 

2000 

C5H5N 

00 

69.99 

18.70 

30.85 

46.56 

53. 14 

56.27 

C5H4N- 

91.7b 

7  2.6 

18.0 

28.9 

42.8 

48.5 

51.1 

ic 

114±3 

79 

23.6 

32.6 

43.1 

48.2 

50.8 

Iic 

142±6 

tll^ 

152*6 

HCCCN 

84. 6** 

56.54 

15.48 

18.66 

22.38 

24.24 

25.19 

HC-CHCN 

97±3 

64.7 

15.9 

19.8 

25.3 

27.8 

29.3 

H2CCHCN 

43.9® 

65.32 

15.26 

20.93 

28.90 

32.66 

34.53 

CVA^ 

99.5 

76.92 

22.53 

30.17 

39.55 

43.82 

45.99 

HC5N 

142 

73.91 

22.45 

27.03 

32.54 

35.21 

36.52 

C-C4H3NCH38 

46.4 

77.5 

20.3 

32.3 

48.6 

58.2 

62.7 

Groups 

Cd-(H)(CN) 

37.6 

37.7 

10.12 

13.47 

17.67 

19.48 

20.38 

Ct-(CN) 

57.7 

31.8 

10.20 

12.11 

14.31 

15.28 

15.71 

76  *  5 

a.  Ref.  33 

b.  For  0- 

and  p-pyridyls, 

add  -1.5  kcal 

mol~* . 

(See  Ref. 

10) 

c.  For  structures  of  these  open-chain  radicals  from  pyridyls,  see  Fig.  2.  A®Hf  ^oo  values 
have  been  estimated  for  all  other  possible  open-chain  radicals  from  pyridyls  and  are  all 
>150*6  kcal  mol”*- 

d.  Ref.  34 

e .  Ref.  35 

f.  (TVA  ■  cyanovlnylacetylene ,  HC  =  C  -  CH  ■  CH  -  CN 


g.  C-C4H3NCH3 


Table  II:  Reaction  Model  for  Pyridine  Pyrolysis 

FORWARD  REACTION  REVERSE  REACTION 


REACTIONS 

loq  A 

n 

E 

loq  A 

n 

E 

REF 

1 

C5H5N=o-C5H4N+H 

15.90 

0.0 

98.0 

14.10 

0.0 

-9.8 

PW 

2 

C5H5N=C5H4N+H 

15.48 

0.0 

102.0 

13.67 

0.0 

-7.3 

PW 

3 

o-C5H4N=C4H4CN 

14.00 

0.0 

25.0 

11.74 

0.0 

-0.3 

est 

4 

C4H4CN=C2H2+HCCHCN 

13.18 

0.0 

37.0 

11.49 

0.0 

3.7 

PW 

5 

HCCHCN=HCCCN+H 

15.90 

0.0 

40.0 

16.14 

0.0 

0.2 

PW 

6 

C5H5N+H=o-C5H4N+H2 

13.48 

0.0 

5.0 

12.29 

0.0 

1.6 

PW 

7 

C5H5N+H=C5H4N+H2 

13.70 

0.0 

9.0 

12.51 

0.0 

4.1 

PW 

8 

C4H4CN=CVA+H 

12.00 

0.0 

40.0 

11.30 

0.0 

3.2 

est 

9 

C5H4N=HCN+C4H3 

14.48 

0.0 

57.0 

9.65 

0.0 

-4.8 

PW 

10 

C4H4CN=HCN+C4H3 

13.70 

0.0 

45.0 

11.13 

0.0 

7.1 

PW 

11 

C4H4CN>=C2H3+HCCCN 

12.70 

0.0 

37.0 

10.98 

0.0 

4.3 

PW 

12 

H+HCCHCN=H2CCHCN 

13.30 

0.0 

0.0 

14.60 

0.0 

104.6 

est 

13 

C2H2+HCN=H2CCHCN 

13.48 

0.0 

40.0 

15.08 

0.0 

80.0 

est 

14 

H+H2CCHCN=HCCHCN+H2 

13.70 

0.0 

8.0 

13.02 

0.0 

7.8 

est 

15 

HCCHCN«C2H2+CN 

17.90 

0.0 

58.0 

16.50 

0.0 

0.4 

PW 

16 

C4H4CN«VA+CN 

14.70 

0.0 

59.0 

13.01 

0.0 

4.6 

PW 

17 

CN+H2=HCN+H 

13.88 

0.0 

0.0 

14.35 

0.0 

17.9 

36 

18 

CN+C2H2=C2H+HCN 

12.70 

0.0 

3.0 

12.45 

0.0 

0.6 

est 

19 

H+H2CCHCN=HCN+C2H3 

13.00 

0.0 

4.0 

11.12 

0:0 

4.4 

est 

20 

H2CCHCN=HCCCN+H2 

14.00 

0.0 

80.0 

13.56 

0.0 

40.0 

est 

21 

CVA+H»=C2H2+HCCHCN 

13.00 

0.0 

4.0 

12.02 

0.0 

7.6 

est 

22 

C4H3+C5H5N=o-C5H4N+VA 

12.30 

0.0 

8.0 

12.47 

0.0 

6.0 

est 

23 

C4H3+C5H5N«=C5H4N+VA 

12.30 

0.0 

10.0 

12.47 

0.0 

6.5 

est 

24 

CH3+C5H5N=o-C5H4N+CH4 

11.70 

0.0 

10.0 

11.98 

0.0 

7.2 

est 

25 

CH3+C5H5N=C5H4N+CH4 

11.70 

0.0 

11.5 

11.98 

0.0 

7.2 

est 

26 

CVA+H*HCN+C4H3 

13.00 

0.0 

4.0 

11.13 

0.0 

2.9 

est 

27 

C2H+HCN=HCCCN+H 

12.30 

0.0 

0.0 

14.19 

0.0 

20.2 

est 

28 

C2H+HCCCN-C4H2+CN 

12.85 

0.0 

3.0 

12.34 

0.0 

0.3 

PW 

29 

CH3CN+H-HCN+CH3 

13.70 

0.0 

2.0 

12.22 

0.0 

10.2 

PW 

30 

H4-C5H5N-C-C4H3NME 

12.30 

0.0 

2.0 

12.12 

0.0 

41.1 

PW 

31 

C-C4H3NME-CH3CN+C3H3 

14.70 

0.0 

60.0 

11.02 

0.0 

7.4 

PW 

32 

C3H3+H*C3H4P 

12.70 

0.0 

0.0 

15.24 

0.0 

89.1 

44 

33 

2C3H3-C2H2+VA 

12.70 

0.0 

0.0 

13.68 

0.0 

38.2 

est 

34 

CH3+C2H-C3H4P 

12.70 

0.0 

0.0 

16.28 

0.0 

115.6 

est 

35 

C4H3*C2H2+C2H 

14.48 

0.0 

57.0 

13.71 

0.0 

1.7 

18 

Table  IK  cont.  ) 


FORWARD  REACTION  REVERSE  REACTION 


REACTIONS 

log  A 

n 

E 

log  A 

n 

E 

REF 

36 

C4H3=C4H2-t-H 

12.00 

0.0 

46.0 

12.37 

0.0 

5.8 

PW 

37 

C2H+H2=H+C2H2 

12.85 

0.0 

0.0 

13.57 

0.0 

20.2 

37 

38 

C2H+C2H2=C4H2+H 

12.00 

0.0 

0.0 

13.14 

0.0 

15.1 

38* 

39 

C2H+C4H2=C6H2+H 

12.70 

0.0 

0.0 

14.03 

0.0 

14.8 

PW 

40 

H+C2H2=C2H3 

12.64 

0.0 

2.4 

12.37 

0.0 

42.8 

39^* 

41 

H-*-C2H4=C2H3+H2 

13.70 

0.0 

8.0 

11.87 

0.0 

6.9 

16' 

42 

CH3+H2=CH4+H 

2.81 

3.0 

7.7 

15.14 

0.0 

17.3 

40 

43 

C2H+C2H3=VA 

13.00 

0.0 

0.0 

16.01 

0.0 

120.8 

est 

44 

C2H4=C2H2+H2 

17.41 

0.0 

79.3 

15.86 

0.0 

37.8 

40 

45 

C4H3+C2H2=L-C6H5 

12.00 

0.0 

3.0 

14.15 

0.0 

36.9 

18 

46 

L-C6H5*C6H5 

10.30 

0.0 

1.4 

13.65 

0.0 

65.2 

18 

47 

C6H6=H+C6H5 

15.70 

0.0 

107.9 

13.10 

0.0 

-2.8 

41 

48 

C2H3+C4H2=L-C6H5 

12.00 

0.0 

3.0 

14.06 

0.0 

36.7 

18 

49 

2C6H5=C12H10 

12.48 

0.0 

0.0 

16.47 

0.0 

107.5 

42 

50 

C6H6+H*C6H5+H2 

14.40 

0.0 

16.0 

12.42 

0.0 

9.8 

43 

51 

2C3H3=L-C6H6 

12.70 

0.0 

0.0 

14.95 

0.0 

57.8 

44 

52 

L-C6H6=C6H6 

12.00 

0.0 

33.0 

16.87 

0.0 

116.0 

44 

53 

C4H4CN+C2H2=C6H6CN 

12.70 

0.0 

1.2 

14.57 

0.0 

37.2 

est 

54 

C6H6CN«C-C6H6CN 

11.30 

0.0 

0.0 

14.22 

0.0 

49.7 

est 

55 

BZN+H-C-C6H6CN 

13.60 

0.0 

2.4 

13.56 

0.0 

26.0 

est 

56 

2C2H2=VA 

13.77 

0.0 

44.6 

15.17 

0.0 

81.1 

26 

57 

VA=H2+C4H2 

14.32 

0.0 

87.0 

13.34 

0.0 

45.4 

26 

58 

C2H+HCCCN-HC5N+H 

12.00 

0.0 

0.0 

11.98 

0.0 

15.3 

est 

Notes ;  Units  for  A  :  cm^,  mol,  s. 

Units  for  E  :  kcal  mol'^ 

PW  indicates  rate  constant  evaluated  in  present  work, 
est  indicates  rate  constant  was  estimated  in  present  work. 

•  Reverse  rate  constant  agrees  with  value  given  in  Reference  at 
1600K. 

Fall-off  value.  Pressure  dependence  as  given  in  Reference. 

®  Revised  vinyl  thermochemistry  presented  in  Reference  has  been  used 
in  the  present  work. 

Species  Identification 

0-C5H4N;  o-pyrldyl,  C5H4N:  m-  and  p-pyridyls,  VA:  vinylacetylene, 
eVA:  cyanovinylacetylene,  BZN:  benzonitrile,  C3H4;  allene  and 
propyne,  C4H4CN:  HC-CH-CH=CHCN. 

"L-"  denotes  open-chain  radical. 

C-C4H3NME; 


O 


Table  III.  Sensitivity  Coefficients  in  the  Pyridine  Systeir  for  Selected  Major  Products 
Evaluated  at  ICOOK  and  at  100  us • 


Reaction  Normalized  Sensitivity  Coefficients 


No . 

C5H5N 

C2H2 

H2 

HCK 

HC3N’ 

C4H2 

C4H4 

CH4 

C6«6 

1 

-0.054 

0.  166 

0.165 

0.165 

0.173 

0. 160 

0.156 

0.393 

0.274 

4 

-0.040 

0.164 

0.125 

0.067 

0.213 

0.290 

0.180 

5 

0. 108 

-0.111 

0.160 

6 

-0.067 

0.254 

0.345 

0.491 

-0.555 

-0.313 

7 

-0.044 

0.117 

0.109 

0.350 

-0. 104 

0.432 

0.229 

-0.277 

9 

-0.015 

0.044 

0.060 

0.063 

10 

-O.OlO 

0.150 

-0.113 

0.191 

0.113 

11 

-0.014 

0.087 

0.140 

0.115 

15 

0.048 

-0.090 

0.111 

-0. 133 

16 

0.225 

22 

0.324 

23 

0.179 

24 

0.175 

28 

0.047 

-0.055 

0.104 

29 

0.734 

30 

0.011 

-0.051 

-0.081 

-0.050 

0.326 

0.351 

31 

0.012 

-0.055 

-0.088 

-0.054 

-0.058 

0.344 

0.378 

33 

0.116 

-0.358 

35 

-0.066 

36 

-0.081 

0.084 

0.354 

-0.185 

37 

0.084 

-0.161 

-0.088 

38 

0.090 

39 

-0.137 

42 

0.260 

45 

0.003 

0.151 

47 

51 

52 
56 


-0.367 


0.078 

0.392 

0.132 


List  of  Figures 


Ti  s  t  r  i  but  ior.  of  all  products  of  significance  behind  reflected 
shocks  in  the  0.7*'  pyridine  in  argon  fixture  over  the 
terperature  range  studied. 
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Abstract 

The  kinetics  of  pyrolysis  of  pyrrole  dilute  in  argon  have  been  studied  in  a  single  pulse  shock  tube, 
using  capillary  column  GC,  together  with  GC/MS  and  FTIR  for  product  identification,  over  the  temperature 
range  1200-1700  K,  total  pressures  of  7.5-13.5  atm  and  nomimal  mixture  compositions  of  pyrrole  of  5000  and 
700  ppm  .  Time-resolved  measurements  of  the  rate  of  disappearance  of  pyrrole  behind  reflected  shock  waves 
have  been  made  by  absorption  spectroscopy  at  230  nm,  at  pressures  of  20  atm  and  mixture  compositions  between 
1000-2000  ppm  pyrrole.  At  the  lower  end  of  the  studied  temperature  range,  the  isomers  of  pyrrole  -  allyl 
cyanide  and  cis-  and  trans-crotononitrile  -  were  the  principal  products,  together  with  hydrogen  qanide  and 
propyne/allene.  At  elevated  temperatures,  acetylene,  acetonitrile,  cyanoacetylene  and  hydrogen  became 
important  products. 

The  rate  of  overall  disappearance  of  pyrrole,  as  measured  by  absorption  spectrometry,  was  found  to 
be  first  order  in  pyrrole  concentration,  with  a  rate  constant  kji,(pyrrole)sl0^^-^*°-’ 
exp(-74.h3.0  kcal  mol'VRT)  s*^  between  1350-1600  K  and  at  a  pressure  of  20  atm.  First  order  dependence 
of  pyrrole  decomposition  and  major  product  formation  was  also  observed  in  the  single  pulse  experiments  over 
the  range  of  mixture  compositions  studied. 

A  74-step  reaction  model  is  presented  and  shown  to  substantially  fit  the  observed  temperature 
profiles  of  the  major  product  species  and  the  reactant  profile.  In  the  model  the  initiation  reaction  is 
postulated  to  be  the  reversible  formation  of  pyrrolenine,  (2H-pyrrole).  Pyrrolenine  can  undergo  ring 
scission  at  the  C2-N  bond  forming  a  biradical  which  can  rearrange  to  form  allyl  cyanide  and  croti^itrile 
or  undergo  decomposition  to  form  HCN  and  CjH^  or  acetylene  and  a  precursor  of  acetonitrile.  The  model 
predicts  an  overall  rate  of  disappearance  of  pyrrole  in  agreement  with  the  experimental  measurements. 
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Introduction 

There  is  considerable  recent  interest  in  the  reaction  pathways  whereby  fuel-bound  nitrogen  evolves 
as  NO,  in  combustion  of  coals  and  other  heavy  fuels  (1,2].  Coals  and  coal-derived  liquids  contain  nitrogen 
in  heteroqclic  structures  containing  pyridine  and  pyrrole  rings  [3,4].  In  furthering  our  understanding 
of  NO,-foraing  reactions  in  combustion  of  heavy  fuels  it  is  important  to  identify  the  volatile  nitrogen- 
containing  precursors  \diich  form  upon  initial  thermal  decomposition  of  the  original  naturally  occurring 
fuel-nitrogen  structures. 

Recently,  we  investigated  [5]  the  products  of  pyrolysis  of  pyridine  and  the  kinetics  of  evolution 
of  nitrogen  products  from  this  molecule  idiich  might  be  considered  as  representative  of  aromatic  nitrogen 
functional  groups  in  coals.  In  addition  to  hydrogen  qanide,  organic  nitriles,  especially  cyanoacetylene, 
were  found  to  be  important  pyrolysis  products  of  pyridine. 

There  have  been  few  previous  studies  of  the  pyrolysis  of  pyrrole  or  its  derivatives,  however.  Apart 
from  some  brief  studies  of  thermal  [6-8]  and  photochemical  [9]  decomposition  of  pyrrole,  the  only  detailed 
kinetic  study  of  the  pyrolysis  of  pyrrole  is  the  recent  work  of  Lifshit:  et  al  [10]  published  during  the 
course  of  our  present  work.  Lifshitz  et  al  studied  the  thermal  decomposition  of  U  mixtures  of  pyrrole  in 
argon  in  a  single  pulse  shock  tube  (SPST)  over  the  temperature  range  of  1050-1450  K  and  at  pressures  of 
about  3  atm.  From  product  distributions  measured  by  the  SPST  technique  they  found  that  qanoalkene  isomers 
of  pyrrole  (especially  allyl  cyanide  and  cis-  and  trans-croto»nitrilel  were  important  products  of 
decomposition.  Lifshitz  et  al  also  deteimined  from  product  analyses  apparent  Arrhenius  parameters  for 
disappearance  of  pyrrole  and  for  formation  of  major  products.  Although  a  mechanism  was  suggested  for  the 
pyrolysis  of  pyrrole,  they  did  not  carry  out  detailed  kinetic  modeling. 

Our  present  work,  also  carried  out  by  SPST,  not  only  extends  the  results  of  Lifshitz  et  al  to  higher 
temperatures  (1700  K),  but  also  develops  and  evaluates  a  detailed  chemical  kinetic  reaction  model  for 
pyrrole  pyrolysis.  In  addition,  we  report  a  time-resolved  ultraviolet  absorption  spectroscopic  kinetic 
Masurement  of  the  direct  decomposition  of  pyrrole  in  highly  dilute  mixtures  in  argon. 


Pjrolyses  vere  carried  out  on  two  single  pulse  shock  tubes,  one  of  diaoeter  3.8  ci  at  United 
Technologies  Research  Center  which  was  used  priaarily  for  SPST  kinetic  studies;  the  other  at  the  University 
of  Sydney  (diaio.  7.6  cat)  was  used  for  real-tiie  ultraviolet  absorption  studies  and  for  pitiduct 
identification  in  conjunction  with  a  Finnegan  GC/HS  at  CSIRO.  The  United  Technologies  Research  Center  shock 
tube  and  ancillary  on-line  GC  analytical  systea  have  been  described  previously  [11].  The  Sydney  tube  was 
equipped  with  pairs  of  sapphire  UV-transaitting  windows.  The  absorption  spectra  at  230  na  of  pyrrole  vapour 
dilute  in  argon  behind  the  reflected  shock  was  aeasured  at  a  distance  of  35  at  froi  the  endwall  using  a  0.3 
a  Spez  Doubleaate  spectroaeter  with  2  na  slits  and  a  150  V  extended  UV  high  pressure  xenon  arc  la^>  as 
source.  The  tiae  constant  of  the  optical  systea,  photoaultiplier,  oscilloscope  and  associated  circuitry 
was  s5  ps. 

Pyrrole  saaples  were  purchased  froa  Aldrich  (stated  purity  99t)  and  froa  Merck  (purity  981). 
Several  bulb-to-bulb  purifications  were  used  to  produce  saaples  of  puriiy  exceeding  99.81  (by  GC) .  The  only 
detectable  iapurity  peak  was  an  unidentified  species  heavier  than  pyrrole  and  not  present  in  the  reaction 
products.  No  difference  between  the  two  sources  of  pyrrole  was  observed. 

For  the  SPST  kinetic  aeasureaents,  aixtures  of  noainal  co^>ositions  of  5000  and  700  ppi  pyrrole  in 
argon  were  prepared  and  stored  in  glass  tanks.  Pressures  and  teaperatures  behind  the  reflected  shock  were 
coaputed  froa  the  aeasured  incident  shock  velocity.  Residence  tiaes  and  quenching  rates  by  the  rarefaction 
wave  were  aeasured  froa  pressure  profiles  recorded  by  Kistler  gauges.  Reflected  shocked  gas  teq>eratures 
ranged  froa  1200-1700  K  and  total  pressures  froa  7.5-13.5  ata.  Residence  tiaes  at  unifon  tei^ratures 
behind  the  reflected  shock  front  ranged  froa  450-750  ps.  Product  analyses  of  the  SPST  runs  on  the  United 
Technologies  Research  Center  tube  were  aade  on  a  Chroapack  CP  Sil  5  c^illary  coluan. 

The  tiM-resolved  UV  absorption  aeasureaents  of  pyrrole  decoq)osition  were  carried  out  on  aixtures 
of  1000-2000  ppa  pyrrole  in  argon  at  teaperatures  between  130C'1700  A,  pressures  of  20  ata  and  residence 
tiaes  up  to  800  ps.  In  the  absorption  studies  teaperatures  were  calculated  froa  aeasureaents  of  both 
incident  and  reflected  shock  velocities. 

Product  identification  and  calibration  were  generally  aade  by  direct  ccaiparison  with  standards. 
The  CjHjCK  isoaeric  products  were  identified  by  GC/MS  and  retention  tiaes  confiraed  with  coaaercial  saaples. 
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Cyanoacetylene  was  identified  both  by  FTIR  spectroscopy  and  GC/MS.  Both  hydrocaitons  and  nitrogen- 
containing  coopounds  were  analysed  by  FID.  Hydrogen  was  detected  by  TCD.  As  reported  previously  [10]  the 
capillary  column  and  sampling  system  eihibited  a  'memory  effect'  for  the  nitrogen-containing  compounds, 
especially  pyrrole  and  its  qaiioalkene  isomers,  as  well  as  acetonitrile  and  acrylonitrile.  Repeatability 
of  injections  of  these  compounds  was  only  *  lOt  Stock  mixtures  of  pyrrole  in  argon  slowly  decreased  in 
concentration  with  storage.  It  was  found  necessary  to  determine  the  concentration  of  pyrrole  in  the  pre¬ 
shock  mixture  just  prior  to  firing  each  shock.  This  was  effected  by  injecting  a  very  small  volume  of  the 
pre-shock  mixture  via  a  sampling  valve  into  the  GC.  Actual  initial  concentrations  of  pyrrole  in  argon 
ranged  from  3700-5200  ppm  for  the  nominally  5000  ppm  mixture  and  from  550-870  ppm  for  the  nominally  700  ppm 
mixture. 


Results 

Product  distributions 

Yields  of  all  products  of  significance  (scaled  by  the  initial  concentration  of  pyrrole)  are  shown 
in  Figure  1  as  a  function  of  temperature  for  a  mixture  of  5000  ppm  pyrrole  in  argon.  The  variation  of 
pyrrole  concentration  with  temperature  is  also  given  in  this  Figure.  There  is  little  difference  in  the 
normalized  product  distributions  between  the  (nominally)  5000  and  700  ppm  mixtures. 

At  U.W  lowest  temperatures  at  tdiich  decomposition  can  be  detected,  the  cyanoalkene  isomers  of 
pyrrole  -  allyl  qanide  and  cis-  and  trans-crotononitrile  -  are  the  major  products.  HCH  is  also  formed  at 
low  temperatures  in  approximately  equal  concentration  with  each  isomer.  Lesser  amounts  of 
proppe/allene  (the  capillary  column  did  not  fully  resolve  these  two  isomers)  are  observed.  With  increase 
in  temperature,  acetylene  becomes  an  increasingly  important  product.  Hydrogen,  too,  is  a  major  product  at 
higher  temperatures.  The  apparent  activation  energies  of  formation  of  the  last  two  products  are  higher  than 
those  of  the  cyanoalkenes  and  HCN  and  C3H4. 

At  elevated  temperatures,  yields  of  cyanoacetylene,  acetonitrile  and  acrylonitrile  rise  steeply.  Other 
products  of  significance  which  have  large  activation  energies  of  formation  are  methane  and  ethylene.  Much 
smaller  amounts  of  ethyl  qanide  and  methylmethacrylate,  the  branched  chain  isomer  of  CjHsCN,  were  observed. 
Small  amounts  of  benzene  and  benzonitrile  are  found  at  high  temperatures.  Traces  of  propene,  1,3-butadiene 
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and  ethane  were  also  detected.  K  very  stall  peak  corresponding  to  the  fonula  04113(1  was  obsened  in  the 
GC/MS.  This  was  attributed  either  to  CH3CCCH  or  to  CH2=C=CHCN.  Diacetylene  and  triacetylene  only  becaie 
itportant  at  teiperatures  above  1600  K.  At  high  tetperatures  a  loss  of  caiton  vas  observed.  This  loss 
increased  vith  increase  in  pyrolysis  teiperature  and  cotprised  about  71  of  the  carbon  in  the  reactant  at 
1700  K.  Sote  tarry  taterial  vas  observed  in  the  high  tei^rature  pyrolysate. 

The  observed  distribution  of  products  is  quite  siiilar  to  those  observed  by  Lifshitz  et  al  [10] 
except  that  the  teiperatures  of  onset  of  products  were  lover  in  the  earlier  studies  [10]  because  of  the 
longer  residence  tiies  (:2  is).  Snail  differences  in  product  distributions  betveen  the  tvo  studies  involve 
the  Minor  product,  ethane  (Lifshitz  et  al  observe  substantially  lore  than  ve)  and  in  the  cis/trans  ratio 
of  crotononitrile.  We  found  this  ratio  to  vary  little  vith  tenperaturc,  ranging  fron  approxiiately  1.5  at 
1300  K  to  1.2  at  1600  K  at  irtiich  temperature  appreciable  deconposition  of  crotononitrile  took  place.  On 
the  other  hand,  Lifshitz  et  al  reported  that  the  product  vas  nearly  all  cis  at  the  lovest  temperatures  of 
their  study  (*1050  K)  vith  the  cis/trans  ratio  decreasing  to  about  1.4  around  1300  K.  Presunably,  the 
cis/trans  ratio  observed  in  the  present  vork  vas  a  consequence  of  rapid  cis-trans  equilibration  taking  place 
at  elevated  teiperatures. 

The  higher  yields  of  ethane  recorded  by  Lifshitz  et  al  probably  arose  froi  enhanced  secondary 
reactions  taking  place  at  their  higher  mixture  concentrations  (11  pyrrole  in  argon) . 

Time-resolved  kinetic  measurements 

An  oscilloscope  trace  of  the  230  ni  absorption  by  a  shock  heated  pyrrole  mixture  (1000  ppm)  is  shovn 
as  a  function  of  time  in  Figure  2.  The  vavelength  studied  corresponds  to  the  lovest  transition  of 
pyrrole  [12,13].  In  Figure  2,  it  may  be  seen  that  there  is  an  increase  in  absorption  due  to  the  increase 
in  density  after  passage  of  the  incident  shock  front.  A  further  large  increase  in  absorption  takes  place 
initially  upon  passage  of  the  reflected  shock,  folloved  by  exponential  decay  of  absorption  resulting  from 
decomposition  of  the  pyrrole  behind  the  reflected  shock  front.  At  the  lov  densities  of  pyrrole  employed 
in  the  absorption  studies,  the  decrease  in  temperature  due  to  reaction  is  less  than  1  X. 

A  plot  of  ln(-ln(I/I„)}  vhich  is  proportional  to  ln{  (pyrrole]  J,  is  shovn  as  a  function  of  time  for 
this  run  in  Figure  3.  (I„  and  I  are  the  unattenuated  intensity  and  intensity  at  tine  t,  respectively.) 
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These  plots  are  linear  except  for  a  snail  curvature  at  late  tines.  First  order  rate  constants  have  been 
obtained  fron  the  initial  linear  portions  of  these  plots.  At  high  tenperatures  (>1650  Kj  the  oscillograns 
of  light  transnission  asynptote  practically  to  I,  indicating  there  are  only  lo¥  concentrations  of  absorbers 
of  230  nn  radiation  in  the  high  tenperature  products.  At  tenperatures  of  1700  K  and  above,  the  rate  of 
decay  of  absorption  at  230  nn  becones  conparable  with  the  tine  constant  of  the  detection  systen. 

An  activation  energy  plot  of  the  first  order  rate  constants  derived  fron  the  absorption  of  230  nn 
radiation  by  pyrrole  is  shown  in  Figure  4.  Regression  analysis  yields  Arrhenius  paraneters  for  the  rate 
constant,  k^^,,  for  overall  disappearance  of  pyrrole  of  log(A/s'M=14.h0.6  and  E,=74.1t3.0  kcal  nol'^ 
between  1480*1680  K. 

Discussion 

Thenochenical  considerations 

As  shown  earlier  [10]  the  najor  pathways  in  the  pyrolysis  of  pyrrole  involve  isonerization  to  the 
nitriles,  C3H5CII,  together  with  fragnentation  principally  into  HCN,  C3H4,  and  CH3CN.  Based  on 
neasurenents  of  apparent  Arrhenius  paraneters  for  fomation  of  the  najor  products,  Lifshit:  et  al  [10] 
suggested  that  the  nitrile  isoners  and  the  HCN  and  C3H4  arose  fron  a  single  transition  state  which  involved 
a  sinultaneous  uninolecular  C-N  bond  cleavage  together  with  a  sinultaneous  H-transfer  of  the  N  hydrogen  to 
the  C2  position.  A  further  H-aton  transfer  in  the  biradical  thus  foned  would  be  necessary  for  fomation 

of  C3H5CN. 

In  the  case  of  furan,  the  0-analogue  of  pyrrole,  pyrolysis  is  initiated  by  C-0  fission  followed  by 
subsequent  rearrangenent  and/or  fragnentation  of  the  biradical  [14-16].  For  furan,  the  rate  constant  for 
disappearance  is  [15,16]  kj^.lfuran)  =  10^®-^eip(-78  kcal  nol'VNT)  s■^  as  expected  for  ring  opening.  In 
contrast,  cur  value  for  pyrrole,  obtained  fron  the  tine-resolved  absorption  neasurenents,  of  kj^,(pyrrole) 
.  ioi<-i*°-6exp(-74.h3.0  kcal  nol'VRT)  s'\  does  not  appear  to  be  representative  of  typical  Arrhenius 
paraneters  for  direct  ring  opening  to  a  biradical. 

A  consideration  of  the  theraochenistry  of  the  pyrrole  systen  can  be  of  assistance  in  unravelling 
the  pyrolysis  nechanisn.  It  is  generally  agreed  that  pyrrole  is  considerably  nore  aronatic  in  character 
than  furan  [17].  Thus,  the  C-N  bond  in  pyrrole  should  be  significantly  stronger  than  a  typical  C-N  single 
bond.  There  are  no  experinental  data  available  for  the  C-N  bond  dissociation  energy  in  pyrrole.  Ue  are 
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forced,  therefore,  to  nake  approxiaate  group  additivity  estiaates  of  the  theraocheaistry  of  the  pyrrole 
systea.  Biradical  (I)  foiaed  on  C-N  scission  vould  have  resonance  foras: 

H  H  H  H  H 

\  /  \  /  \ 

c - C  C - C  C- 

//  W  ^  /  •  ^  / 

M  ■  N  ^ 

V  \  \ 

H  H  H 

(D  (ID  (IID 

An  approxiaate  estiaate  of  iH^.joofn)  =  123+5  kcal  aol'^  can  be  obtained  froi  groups  Hj-IH),  Cj-(H)(Cj), 
C  ,j-(H)(Cj)  and  C  -(H)(Cj)2,  available  in  the  literature  [18,19],  The  biradical  will  be  resonance 
stabilized.  A  stabilization  energy  of  about  7  kcal  aol'\  siailar  to  that  in  the  biradical  foned  during 
furan  decoaposition  [14,16],  would  lead  to  an  enthalpy  of  fonation  of  the  stabilized  biradical  of  about 
116i6  kcal  aol*^  and  to  a  C-H  bond  dissociation  energy  in  pyrrole  of  approxiaately  90  kcal  aol*^  If  these 
estiaates  are  correct,  then  our  obsened  activation  energy  for  disappearance  of  pyrrole  of  74  kcal  aol‘\ 
or,  indeed,  that  reported  by  Lifshitz  et  al  (75  kcal  aol*^)  would  not  appear  to  correspond  to  direct  ring 
scission. 

Pyrroles  are  known  to  rearrange  theiaally  at  teaperatures  between  500-60(K  [20]. 

Patterson  et  al  [21-23]  studied  the  theraal  rearrangeaent  of  several  A-substituted  pyrroles  to  the 
corresponding  2-substituted  pyrrole.  Arrhenius  paraaeters  for  several  H-alkylpyrrole  rearrangeaents  have 
been  aeasured  by  Jacobson  and  Jensen  [24,25].  The  rearrangeaents  were  found  to  be  hoaogeneous 
intraaolecular  shifts.  Pyrrolenine  (2H-pyrrole)  interaediates  were  postulated  to  be  involved  in  the 
aechanisa.  Pyrrole  itself  undergoes  this  rearrangeaent.  Patterson  et  al  [23]  pyrolysed  N-deuterated 
pyrrole  under  siailar  conditions  to  the  B-alkyl  rearrangeaents  and  found  2-deuterated  pyrrole.  Although 
pyrrolenines  have  not  been  isolated  froa  the  pyrolysis  products  of  siaple  K-substituted  pyrroles,  2H- 
pyrroles  have  been  obsened  in  the  photolysis  of  H-alkyl  ring  substituted  pyrroles  [20].  There  is  now 
considerable  indirect  evidence  for  the  presence  of  pyrrolenine  interaediates  in  both  theraal  and 
photocheaical  rearrangeaents  of  IH-pyrroles. 
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Unlike  its  iN-isooer,  pjrrolenine  would  be  expected  to  possess  a  relatively  weak  C-N  bond  with  the 
C2  carbon.  Again,  there  are  no  eiperinental  themocheaical  data  available  on  pyrrolenine  and  we  have  been 
forced  to  resort  to  approiioate  group  additivity  nethods  to  estinate  the  enthalpy  of  foraation  of  this 
Bolecule.  Under  the  assunption  that  the  ring  correction  tens  are  the  sane  as  for  cyclopentadiene,  we 
estinate  that  4H,_3oo=45.5  kcal  nol'^  for  pyrrolenine.  This  is  in  satisfactory  agreenent  with  the  value  of 
44.6  kcal  nol'^  obtained  by  a  HINDO/3  calculation  of  the  pyrrole/pyrrolenine  energy  of  tautonerization  [26]. 

We  suggest  that  the  pyrolysis  of  pyrrole  is  initiated  by  the  reversible  hydrogen  shift  to 
pyrrolenine: 


(1) 


Pyrrolenine  nay  then  undergo  themolysis.  It  would  be  expected  that  the  ring  would  open  by  scission  of  the 
C2-N  bond. 


H  H 

\  / 
C=C  H 

(IV) 


Biradical  (IV)  could  then  undergo  rapid  H  aton  shifts  to  fora  cis-crotononitrile  or  allyl  cyanide.  The 
first  isoner  can  be  foraed  by  a  sinple  1,4  H*shift  although  to  fora  the  second  would  appear  to  require  a 
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1,2  H-shift,  usually  a  process  of  relailvely  high  activation  energy.  However,  as  pointed  out  by  Benson 


[27],  this  shift  can  be  facile  and  rapid  in  a  biradical. 

Hydrogen  qanide  and  propyne/allene  could  arise  fros  a  transition  state  resulting  froi  contraction 
of  the  biradical  to  a  cyclopropenyl  ring  with  subsequent  eliiination  of  HCN.  Propyne  and  allene  would  then 
be  the  products  of  isoaerization  of  a  'hot'  cyclopropane  lolecule  [28]. 

In  a  process  of  higher  activation  energy,  the  biradical,  having  undergone  1,3-hydrogen  shifi,  could  then 
elininate  acetylene,  foning  a  C2H3R  precursor  of  CH3CR  such  as  2H-aziridine  (V)  which  would  then  rapidly 
rearrange  to  acetonitrile.  This  rearrangeaent  has  been  studied  theoretically  by  ab  initio  aethods  [29,30] 


CH3CN 

To  consider  the  feasibility  of  this  suggested  aechanisa  it  is  necessary  to  have  an  estiaate  of  the 
enthalpy  of  foraation  of  biradical  (IV).  This  biradical  should  possess  allylic  resonance  stabilization. 
An  approxiaate  value  of  the  enthalpy  of  fonation  can  be  obtained  using  groups  froa  the  allyl  radical 
together  with  groups  Cj-(H)(Cj),  Cj-(H)(!(J  and  Tj.  The  last-aentioned  group  contribution  is 

estiaated  [5]  by  coaparing  group  differences  between  H2C=C*H  and  H2C=CH2  with  H2C=r  and  HjC^HH.  In  this 
way  we  estiaate  AH,  30^1 IV)  =  118+5  kcal  aol'^ 

Thenocheaical  paraaeters  used  in  the  present  study  and  not  available  in  the  theraocheaical 
literature  are  presented  in  Table  I. 
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Reaction  model 

The  detailed  reaction  Dodel  to  simulate  the  fields  of  the  major  products  is  given  in  Table  II.  As 
suggested  above,  initiation  is  via  the  reversible  reaction  (1)  to  form  pyrrolenine.  Although  activation 
energies  for  N-alkyl  rearrangements  in  pyrroles  are  between  57+3  kcal  mol*^  and  A-factors  around  10“s'^ 
[23,25],  no  kinetic  data  are  available  for  H-H  to  C2  shifts.  We  have  therefore  chosen  to  optimize  the  rate 
parameters  of  reaction  (1).  As  discussed  below,  the  reaction  model  is  very  sensitive  to  the  value  of  k^. 
Optimization  involved  the  systematic  variation  of  A^  and  to  best  fit  the  experimental  pyrrole  profile 
and  profiles  of  the  major  products.  The  optimized  value  of  ki=3xl0“exp(-67  kcal  mol‘^/RT)  s’^  was 
obtained.  Reactions  (2),  (3),  (4)  and  (5)  of  Table  II  are  not  elementary  reactions  but  involve  the 
intermediate  biradical  (IV)  above,  whose  subsequent  decomposition  is  considered  to  be  rapid  by  parallel 
pathways,  one,  a  rearrangement  to  the  CjH^CN  isomers,  the  second,  a  fission  to  HCX  and  proppe/allene 
(probably  via  a  cyclopropene  intermediate),  rfiereas  the  third  involves  the  elimination  of  C^H^. 

The  same  optimized  value  for  the  activation  energies  of  reactions  (3)  and  (4)  was  obtained,  68  kcal 
mol'\  This  activation  energy’  is,  within  the  error  of  estimating  the  enthalpies  of  formation,  equal  to  the 
enthalpy  difference  between  biradical  (IV)  and  pyrrolenine.  For  purposes  of  modeling,  both  cis-  and  trans* 
crotononitrile  have  been  lu^ed  together.  A-factors  for  reactions  (3)  and  (4)  are  5x10“  s*S  normal  values 
for  ring  scission  to  a  biradical  [15,16] .  The  optimized  rate  constant  for  reaction  (2)  is  k2=3xl0“exp(-56 
kcal  mol'VRT)  This  preexponential  factor  is  suggestive  of  a  qclic  transition  state  and  the 
activation  energy  is  sufficiently  high  to  permit  a  qclopropene  intermediate.  In  the  model,  has  been 
considered  to  be  propyne,  the  more  thermodynamically  stable  form  at  high  temperatures.  The  rate  constant 
for  reaction  (5)  derived  from  the  modeling  is  k5=6xl0“exp(-76  kcal  mol’VRT)  s'^  The  derived  activation 
energy  is  probably  high  enough  to  enable  biradical  (IV)  to  split  into  CjHj  and  a  C2H3N  precursor  of 
acetonitrile  such  as  2H-aziridine  although  other  intermediates  of  lower  energy  such  as  ketenimine  have  been 
identified  on  the  surface  [29,30). 

The  model  contains  secondary  reactions  of  the  major  products.  Radicals  important  in  the  mechanism 
can  be  identified  by  their  products  of  termination.  Presumably,  ethane  arises  from  methyl  radical 
recombination.  Ethyl  cyanide  results  from  combination  of  CH3  and  CHjCN.  A  important  reaction  for 
production  of  the  latter  radical  is  the  unimolecular  fission  reaction  of  allyl  cyanide  [reaction  (-32)). 
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HoBogeneous  tensination  of  vinyl  radicals  would  be  expected  to  produce  1,3-butadiene.  The  presence  of  a 
C^HjN  product  (either  CH3CCCN  or  CH2=C=CHC1I)  suggests  the  intenediacy  of  a  cyanoallyl  radical.  This 
radical  can  be  produced  by  H-abstraction  reactions  of  allyl  cyanide. 

The  observed  yields  of  hydrogen  and  lethane  could  not  be  siBulated  Berely  froB  secondary  reactions 
of  product  species.  To  Bodel  the  ezperiBental  and  CH^  it  was  necessary  to  include  abstraction  reactions 
of  pyrrole  in  the  sechanisB.  The  weakest  bond  involving  hydrogen  should  be  the  N-H  bond  which  Bight  be 
expected  to  be  soBev^at  weakened  owing  to  the  partial  aroBaticity  of  the  CNC  bonds.  The  pyrrolyl  radical 
foned  froB  H  and  CH3  abstraction  of  the  N-hydrogen  is  analogous  to  the  cyclopentadienyl  radical  and  should 
possess  considerable  resonance  stabilization.  Hn  estiaate  of  the  enthalpy  of  fonation  of  pyrrolyl  can  be 
Bade  under  the  assuaption  that  the  ring  stabilization  energy  is  siailar  to  that  in  cyclopentadienyl.  The 
value  thus  obtained  is  AH,  3oo(pyrrolyl)  =  54.4  kcal  aol'^  The  H-bond  fission  reaction  of  pyrrole 
[reaction  (56)]  has  been  included  in  the  aodel  together  with  H  and  CH3  abstraction  reactions  (59,60). 

Ring  opening  reactions  of  the  stabilized  pyrrolyl  radicals  are  likely  to  have  high  activation 
energies.  These  radicals,  like  qclopentadienyl  radicals  [32],  should  be  unreactive.  In  the  Bodel, 
teiBination  reactions  of  these  radicals  to  fon  a  bipyrrolyl  product  have  been  included.  Uhilst  this 
product  has  not  been  specifically  identified,  experiientally  there  was  a  loss  of  carbon  and  production  of 
heavy  product  which  sight  have  contained  soae  nitrogen. 

Integration  of  the  kinetic  equations  of  Table  II  was  Bade  using  the  CHEKKIN  package  [33],  LSODE  [34] 
and  the  SAKDIA  shock  tube  code  [35]  aodified  to  account  for  the  effects  of  quenching  by  the  reflected 
rarefaction  wave.  Rate  of  production  and  sensitivity  analyses  were  carried  out  by  the  SENKIN  code  [36]. 

CoBparison  of  the  kinetic  Bodel  with  experisent  is  given  in  Figures  5-11  for  the  Bajor  products. 
The  Bodel  adequately  describes  the  teiperature  dependence  of  these  products  for  both  series  of  runs  (5000 
and  700  ppi).  The  lodel  also  reproduces  the  experiBental  pyrrole  profiles  up  to  about  801  decoaposition 
(Figure  12).  The  residual  pyrrole  concentration  at  high  tesperatures  results  froB  boundary  layer  effects 
in  the  shock  tube.  The  failure  of  the  Bodel  to  sisulate  the  experiBental  tailing  off  in  decosposition  of 
the  C3H,CN  isoBers  and  the  C3H4  at  high  teaperatures  is  also  probably  a  consequence  of  boundary  layer 
effects. 
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A  brief  selection  of  the  results  of  sensitivity  analyses  of  the  reaction  lechanisa  is  given  in  Table 
III.  This  table  contains  nonaalized  sensitivity  coefficients  evaluated  at  1398  R  and  at  100  iis  for  the  5000 
ppo  aixture  in  argon.  The  sensitivity  analysis  shows  that  the  reactions  Bost  inportant  to  the  decoi^osition 
of  pyrrole  and  fonation  of  najor  products  are  reactions  (l)-{5),(32),(56)  and  (59).  The  first  five 
reactions  all  involve  fonation  or  decoiposition  of  pyrrolenine.  Reaction  (32)  which  produces  vinyl 
radicals  fron  allyl  qanide  is  iaportant  in  detenining  the  production  of  H  atons  and  acetylene. 
Decooposition  of  pyrrole  via  H  atoa  fission  (56)  and  H-abstraction  (59)  is  a  relatively  linor  pathway.  The 
Arrhenius  paraaeters  of  kjg  were  obtained  by  optiaization  to  low-teaperature  eiperiaental  profiles  which 
show  sensitivity  to  this  rate  constant.  The  value  thus  obtained  of  1^^=  10^^  ’  s‘^  This  value  aight 
appear  to  be  soaewhat  low  for  a  uniaolecular  H-fission  reaction.  Increase  in  by  a  factor  of  two  to 
10^^-°  would  still  predict  the  low-teaperature  Hj  profiles  and  the  pyrrole  decay  profiles  within 
eiperiaental  error,  whereas  an  increase  in  A^^  to  10^*-^  would  significantly  overpredict  the  low  teaperature 
yields  of  hydrogen  and  rate  of  pyrrole  decoaposition.  However,  the  A-factor  for  reaction  (56)  light  well 
be  abnoraally  low  -  the  transition  state  probably  reseables  the  tight  resonance  stabilized  pyrrolyl  radical. 

It  reiains  to  interpret  the  eiperiaental  rate  constant  k^,,  for  overall  disappearance  of  pyrrole. 
The  reaction  aodel  was  run  for  teaperatures  between  1300-1600  K.  Over  the  entire  unifora  residence  tiae 
the  aodel  predicted  that  the  rate  of  disappearance  of  pyrrole  obeyed  kinetics  which  were  first  order  in 
pyrrole  concentration.  Linear  regression  of  an  Arrhenius  plot  of  the  theoretically  derived  rate  constants 
93’®  l'di..«od.i(Pyrrole)=10^*’’°“P('’^7.2  kcal  aol'VST)  The  aodel  thus  yields  a  rate  constant  for 
overall  disappearance  of  reactant  which  lies  within  the  error  liiits  of  the  rate  constant  aeasured  by  tiae- 
resolved  ansoiption  spectroaetry  at  230  na.  These  values  aay  also  be  coapared  with  an  initiation  rate 
constant  obtained  by  Lifshitz  et  al  [10]  of  kj„jj(pyrrole):10^*  ®^eip(-75  kcal  aol'VRT)  s’^  Ho  error 
liaits  were  stated  by  these  authors.  However,  since  their  values  were  derived  in  the  liait  of  low  extents 
of  decoaposition  tdiere  errors  in  rate  paraaeters  can  result  due  to  low  product  yields,  it  would  appear  that 
our  Arrhenius  paraaeters  agree  with  their  initiation  rate  constant  within  reasonable  eiperiaental  errors. 
Lifshitz  et  al  also  present  a  rate  constant  for  overall  disappearance  of  pyrrole  of  k„„„ii:10“-®eip(-80 
kcal  aol'VRT)  obtained  froa  aeasureaents  of  (pyrrolej^  (t=t„J  at  different  teaperatures.  At  1300  K  their 
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value  of  ko,,r.u  is  a  factor  of  five  larger  than  our  value  of  k^i,.  The  reason  for  this  disagreeaent  is 
uncertain.  Our  value  of  k^^.  has  been  obtained  in  tvo  different  laboratories  and  by  tvo  different 
techniques.  Possibly  the  use  by  Lifshitz  et  al  of  aore  concentrated  pyrrole  aixtures  leads  to  a  greater 
iaportance  of  biaolecular  decoaposition  of  pyrrole.  Another  possibility  is  the  neglect  of  reaction  in  the 
cooling  transient  in  the  analysis  aade  by  Lifshitz  et  al.  (Reaction  through  the  rarefaction  cooling  was 
taken  into  account  in  our  aodeling  studies.) 

Ue  have  also  unsuccessfully  tested  other  pyrolytic  aechanisas.  In  particular,  we  have  been  unable 
to  fit  a  H'Chain  aodel  to  the  ezperiaental  data.  In  that  aodel  rate  constants  for  R-H  fission  [Reaction 
(56)1  and  for  abstraction  fron  pyrrole  by  H  atcas  [Reaction  (59)]  were  increased  to  siiulate  the 
experinental  decay  of  pyrrole.  It  was  possible  to  postulate  a  H-chain  decowposition  of  pyrrole  by  invoking 
the  intemediacy  of  a  hydrogen  adduct  to  pyrrole  which  could  then  rearrange  and  deconpose  to  the  pyrrole 
isoners,  allyl  qanide  and  crotononitrile,  regenerating  K  atODS.  This  nodel  greatly  over-predicted  the 
experinental  yields,  however.  Additionally,  if  kj^  and  k,,  were  adjusted  to  fit  the  experinental  pyrrole 
profiles  for  the  5000  ppn  nixture,  the  nodel  greatly  underestinated  the  experinental  pyrrole  deconposition 
rate  in  the  700  ppn  nixtures.  The  experinental  species  yields  at  the  two  different  nixture  conpositions 
cannot  both  be  sinulated  by  a  binolecular  nechanisn  of  pyrrole  deconposition. 

A  nodel  in  which  the  pyrrole  deconposed  via  hydrogen  abstraction  reactions  with  fomation  of 
pyrrolyl  but  without  hydrogen  addition  to  pyrrole,  was  also  tested.  It  is  necessary  to  postulate  ring 
opening  of  pyrrolyl  followed  by  addition  of  H  atons  and  abstraction  of  H  fron  all  donors  to  fon  the  pyrrole 
isoneric  products,  crotononitrile  and  allyl  qanide.  This  nodel  was  also  incapable  of  predicting  the 
experinental  species  profiles  at  both  nixture  conpositions  studied. 

On  the  other  hand,  our  pyrolysis  nodel  involving  pyrrolenine,  although  sonewhat  speculative,  is  in 
accord  with  the  pyrrole  tnemochenistq  and  is  capable  of  reproducing  the  experinental  product  distributions 
and  rates.  Our  reaction  nodel  has  sinilarities  to  the  pyrolysis  schene  suggested  by  Lifshitz  et  al  [10]. 
Although  they  did  not  invoke  the  intemediate  pyrrolenine,  their  nechanisn  did  involve  the  sane  transition 
state  for  fomation  of  the  CjIijCN  isoners  as  well  as  acetylene  and  CjH^. 


Conclusions 


The  pffolfsis  kinetics  and  distribution  of  najor  products  can  be  satisfactorily  lodeled  by  a 
reaction  nechanisB  in  idiich  pyrolysis  is  initiated  by  the  reversible  fomation  of  pyrrolenine  (2H-pyrrole). 
Pyrrolenine  is  postulated  to  undergo  C2-N  scission  to  fon  a  biradical  idiich  can  rearrange  to  allyl  cyanide 
and  crotononitrile.  The  sane  biradical  can  decospose  into  HCN  and  propyne/allene  (probably  via  a 
qclopropene  intenediate)  or  into  acetylene  and  a  precursor  of  acetonitrile  such  as  ZH-aziridine  or 
keteniffline.  Experiientally  it  was  found  that  the  overall  rate  of  disappearance  of  pyrrole  is  first  order 
in  pyrrole  concentration.  The  rate  constant  for  overall  disappearance  of  pyrrole  derived  froi  the  lodel 
is  in  agreeoent  within  experiiental  error  with  an  experiiental  detenination  by  tiie-resolved  absorption 
spectrometry  and  with  a  value  obtained  earlier  by  Lifshitz  et  al  [10]. 
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Figure  1; 

Figure  2: 


Figure  3: 

Figure  4: 

Figures  5-11: 


Figure  12: 


Captions  to  Figures 

Teiperature  distribution  of  all  products  of  significance  in  the  pyrolysis  of  a  5000  ppi 
fixture  of  pyrrole  in  argon.  Total  pressure  13  ati,  residence  tiie  400  ps. 

Oscillograa  of  the  transaission  of  light  of  wayelength  230  no  in  a  shock  vave  heated 
lixture  of  1000  ppo  pyrrole  in  argon.  Reflected  shocked  gas  te^>erature  =  1513  K,  pressure 
=  20  ato.  ISF:  incident  shock  front,  RSF:  reflected  shock  front.  Positions  of  zero  and 
unattenuated  light  intensity  are  indicated. 

Plot  of  ln(I/I„)  versus  tioe  for  the  run  shown  in  Figure  2.  The  rate  constant  derived  froo 
this  run  was  kdi,(1513  K)  =  2.35Il0^s■^ 

irilienius  plot  of  the  first  order  rate  constant,  kj^.,  for  overall  disappearance  of  pyrrole 
as  leasured  by  tioe  resolved  absorption  spectrooetry  at  230  nn. 

Conparison  of  the  experinental  product  yields,  Y,  of  esignated  species  with  the 
predictions  of  the  reaction  Bodel  of  Table  II.  o  and  0:  data  points  for  fixtures  of  5000 

ppf  pyrrole  in  Ar,  x  and  •:  data  points  for  fixtures  of  700  p{«  pyrrole  in  It. - sodel 

predictions  for  5000  ppf.  •  -  •  fodel  predictions  for  700  ppf. 

Coaparison  between  the  cxperifental  pyrrole  profiles  and  the  predictions  of  the  reaction 

fodel  of  Table  II.  0  and  - :  experifental  points  and  fodel  predictions  for  5000  ppa 

fixture  of  pyrrole  in  Ar;  x  and  -  -  -  experifental  points  and  aodel  predictions  for  700 
ppa  fixture  of  pyrrole  in  Ar. 


Table  I:  TheniocheiDical  Parameters  for  the  Pyrrole  System 
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Table  III:  Sensitivity  Analysis  for  the  Reaction  Model  of  Table  II. 

Entries  are  the  nonalised  lass  fraction  sensitivities  for  the  given  species  evaluated  at  100  ps  and  tenperature  1398  K. 


Table  III  (cont.) 
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